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ABSTRACT 

We present a systematic study of the X-ray and multiwavelength properties of a sample of 17 highly radio- 
loud quasars (HRLQs) at z > 4 with sensitive X-ray coverage from new Chandra and archival Chandra, 
XMM-Newton, and Swift observations. Eight of the new and archival observations are reported in this work 
for the first time. New Chandra observations of two moderately radio-loud and highly optically luminous 
quasars at z > 4 are also reported. Our HRLQ sample represents the top ~ 5% of radio-loud quasars in terms 
of radio loudness. We found that our HRLQs have an X-ray emission enhancement over HRLQs at lower 
redshifts (by a typical factor of sa 3), and this effect, after controlling for several factors which may introduce 
biases, has been solidly estimated to be significant at the level. HRLQs at z — 3-4 are also found to 

have a similar X-ray emission enhancement over z < 3 HRLQs, which supports further the robustness of our 
results. We discuss models for the X-ray enhancement's origin including a fractional contribution from inverse 
Compton scattering of cosmic microwave background photons. No strong correlations are found between the 
relative X-ray brightness and optical/UV emission-line rest-frame equivalent widths (REWs) for radio-loud 
quasars. However, the line REWs are positively correlated with radio loudness, which suggests that relativistic 
jets make a negligible contribution to the optical/UV continua of these HRLQs (contrary to the case where the 
emission lines are diluted by the relativistic ally boosted continuum). Our HRLQs are generally consistent with 
the known anti-correlation between radio loudness and X-ray power-law photon index. We also found that the 
two moderately radio-loud quasars appear to have the hardest X-ray spectra among our objects, suggesting that 
intrinsic X-ray absorption (Nu ~ 10 23 cm~ 2 ) may be present. Our z > 4 HRLQs generally have higher X-ray 
luminosities than those for the composite broad-band spectral energy distributions (SEDs) of HRLQs at lower 
redshift, which further illustrates and supports the X-ray emission enhancement of z > 4 HRLQs. Some of our 
HRLQs also show an excess of mid-infrared emission which may originate from the synchrotron emission of 
the relativistic jets. None of our z > 4 HRLQs is detected by the Fermi LAT two-year survey, which provides 
constraints on jet-emission models. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: high-redshift — jets: galaxies — X-rays: 
galaxies 



1. INTRODUCTION 

X-ray studies of z > 4 active galactic nuclei (AGNs) over 
the past w 12 yr have greatly improved understanding of the 
growth of the first supermassive black holes (SMBHs) in the 
early Universe, thanks to unprecedentedly powerful X-ray ob- 
servatories as well as wide-field surveys in the optical/UV 
and radio bands. Chandra and XMM-Newton, along with 
other X-ray missions, have observed more than 100 quasars 
at z > 4, the majority of which were first identified by wide- 
field optical/UV surveys, e.g., the Sloan Digital Sky Survey 
(SDSS; York et al. 2000), or radio surveys, e.g., the Faint Im- 
ages of the Radio Sky at Twenty-Centimeters (FIRST; Becker 
et al. 1995) survey and the NRAO VLA Sky Survey (NVSS; 
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Condon et al. 1998). The quasar population is divided into 
radio-quiet and radio-loud subcategories based on the appar- 
ently bi-modal distribution of the radio-loudness parameter 
(e.g., Ivezic et al. 2004)0 which was defined by Kellermann 
etal. (1989) asR~fs ghz//44oo A (rest-frame; see ij3]for a de- 
tailed definition). Radio-quiet quasars (RQQs) have R < 10, 
while radio-loud quasars (RLQs) are those with R ^ 10. 

The basic X-ray properties of z > 4 RQQs are now well es- 
tablished (e.g., Kaspi et al. 2000; Brandt et al. 2001, 2002; 
Vignali et al. 2001, 2003a, 2003b, 2005; Shemmer et al. 
2005, 2006). The z > 4 RQQs have remarkably similar X-ray 
properties to those of RQQs at lower redshift, which is no- 
table given the strong evolution of the quasar population over 
cosmological timescales (e.g., Richards et al. 2006b; Croom 
et al. 2009; Jiang et al. 2009). These results suggest the mode 
of SMBH growth via accretion does not strongly evolve with 
redshift. 

RLQs feature powerful jets that originate near the SMBH 

7 There is an ongoing debate regarding the bi-modality of the radio- 
loudness distribution of quasars (e.g., Kellermann et al. 1989; Ivezic et al. 
2002, 2004; Cirasuolo et al. 2003; Singal et al 2011, 2012). Some recent 
studies with new approaches, e.g., investigating the loci of quasars in the op- 
tical/UV spectroscopic parameter space (Zamfir et al. 2008) or the quasar 
radio luminosity function (Kimball et al. 2011b) found two distinct popula- 
tions with different radio properties, supporting the dichotomy of radio-loud 
vs. radio-quiet quasars (but also see, e.g., Broderick & Fender 201 1). 
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and may extend up to hundreds of kpc. These jets possess bulk 
relativistic velocities on at least sub-kpc scales and generate 
strong radio emission that is particularly pronounced in low- 
inclination (highly beamed, foreshortened, core-dominated) 
RLQs. RLQs usually show an enhancement of nuclear X-ray 
emission compared to RQQs, likely due to a contribution from 
the spatially unresolved jet (e.g., Worrall et al. 1987; Miller 
et al. 2011; and references therein). Previous X-ray studies 
of representative RLQs at z>4 (R sa 40-400; e.g., Bassett et 
al. 2004; Lopez et al. 2006; Saez et al. 2011) have shown 
that the factor of X-ray emission enhancement of these ob- 
jects is similar to that of low-redshift RLQs of similar radio 
loudness. These results limit any redshift dependence of the 
nuclear X-ray emission in RLQs (consistent with earlier stud- 
ies to moderate redshifts; e.g., Worrall et al. 1987). The lack 
of strong redshift evolution argues against a dominant X-ray 
generation mechanism for the nuclear jet-linked contribution 
whereby relativistic jet electrons upscatter cosmic microwave 
background photons into the X-ray band (the 1C/CMB model; 
e.g., Tavecchio et al. 2000; Celotti et al. 2001). While 
photon fields associated with the quasar are expected to be 
more relevant on small scales (« 0.1-2 kpc), for z = 4 core- 
dominated RLQs (for which 1" is 7.1 projected kpc and typ- 
ical inclinations of 5-7° correspond to deprojected lengths 
~10 times greater) the spatially unresolved X-ray emission 
potentially includes contributions from the jet to tens of kpc. 
The IC/CMB model is often proposed to explain the X-ray 
knot emission from low-redshift large-scale RLQ jets studied 
with Chandra, but here too the consequent predicted strong 
increase with redshift in the X-ray luminosities of jet knots 
(absolute or relative to jet radio or non-lC/CMB core X-ray 
emission; e.g., Schwartz 2002) is to date not observed (e.g., 
Bassett et al. 2004; Lopez et al. 2006; Marshall et al. 2011; 
Saez et al. 2011). While some degree of IC/CMB emission 
is required, the relative contribution (and associated spatial 
scales) to the observed X-ray emission is currently unclear. 

Some z > 4 quasars with the highest values of radio loud- 
ness (R > 400) have been studied individually in X-rays (e.g., 
PMN J0525-3343, Worsley et al. 2004a; Q 0906+6930, 
Romani 2006; RX J1028.6-0844, Yuan et al. 2005; 
GB 1428+4217, Worsley et al. 2004b; GB 1508+5714, 
Siemiginowska et al. 2003, Yuan et al. 2003, 2006). These 
objects are often referred to as "blazars" due to their ex- 
treme radio loudnesses and other properties^ Most of the 
previous studies were focused on their X-ray spectral prop- 
erties. These objects often show soft X-ray spectral flat- 
tening; the fraction of objects with such flattening appears 
to rise with redshift (e.g., Fiore et al. 1998; Page et al. 
2005; Yuan et al. 2006). This flattening has been pro- 
posed to be caused by intrinsic X-ray absorption with a col- 
umn density of Nh ~ (1-3) x 10 22 cm~ 2 (but also see Be- 
har et al. 2011). This level of intrinsic X-ray absorption 
is also found in one z > 4 quasar with moderate radio loud- 

8 The term "blazar", including subcategories of BL Lac objects and op- 
tically violently variable quasars (e.g., see §1.3 of Krolik 1999), has been 
defined based on a set of criteria, including the equivalent widths of emission 
lines, the non-thermal jet-linked contribution to the optical continuum, broad- 
band variability, the radio spectral shape, and the polarization. However, this 
terminology has sometimes been applied either inconsistently or based on 
incomplete information. For our objects, it is often not feasible to classify 
them clearly as "blazars" based on available data. All of them are broad-line 
quasars (i.e., not BL Lac objects). All but one (J0913+5919) of our objects 
have flat radio spectra (see g3j. In this work, we have made minimal use of 
the term "blazar". Instead we define a highly radio-loud quasar sample based 
on radio loudness (see ij2]for details). 



ness (SDSS J001 1 + 1446; Saez et al. 2011). An alternative 
mechanism for the soft X-ray spectral flattening is the bulk 
Comptonization of broad-line photons by relativistic plasma 
moving along the RLQ jet (e.g., Begelman & Sikora 1987; 
Sikora et al. 1997). Celotti et al. (2007) studied the time- 
dependent spectra generated by this mechanism and repro- 
duced the soft X-ray flattening feature in the XMM-Newton 
spectrum of GB 1428+4217. Volonteri et al. (201 1) studied 
hard-X-ray-selected high-redshift blazars and found a deficit 
of RLQs at z > 3 compared to the expectations from the num- 
ber of blazars. Possible explanations of this deficit require 
either a decrease of average bulk Lorentz factor in relativis- 
tic jets or an increase by a factor of 2-10 in the number of 
very massive black holes (Mbh > 10 9 M@) at high redshifts 
(derived from the radio luminosity function). They stressed 
the importance of finding high-redshift blazars which could 
provide constraints on the SMBH-dark halo connection and 
structure formation. 

Although these highly radio-loud quasars (HRLQs here- 
after) are rare, they produce among the most powerful rela- 
tivistic jets in the Universe. These remarkable objects have 
required the larger sky coverage of modern wide-field opti- 
cal/UV and radio surveys to be discovered in sufficient num- 
bers to assemble a statistically meaningful sample at high 
redshift. Current major X-ray missions provide coverage 
of the 2^4-0 keV energy band in the rest frame of z > 4 
HRLQs. Studies of the X-ray and broad-band properties of 
high-redshift HRLQs provide probes of the jet-launching and 
radiation mechanisms in the early stages of cosmic evolution 
(e.g., Fabian et al. 1998, 1999). However, there has not been a 
systematic survey to investigate the general properties of these 
remarkable HRLQs at z > 4. Miller et al. (201 1; Ml 1 here- 
after) studied the X-ray properties of a large, diverse sample 
of RLQs0 with wide ranges of redshift and radio loudness. 
However, HRLQs at z > 4 only made up < 1 % of their full 
sample. In this work, we conduct the first X-ray survey of a 
sample of HRLQs at z > 4 using new X-ray observations by 
Chandra and sensitive archival X-ray data. Our science goals 
include the following: (1) assess any X-ray emission enhance- 
ment of these objects and compare them to HRLQs at lower 
redshift; (2) provide constraints on the basic X-ray spectral 
properties of the newly observed Chandra targets, and study 
the relation between X-ray spectral properties and radio loud- 
ness; (3) investigate relations between relative X-ray bright- 
ness, radio loudness, and optical/UV emission-line proper- 
ties; and (4) study the broad-band spectral energy distribu- 
tions (SEDs) of our objects. 

In $2] we describe the selection of our sample of HRLQs 
at z > 4. In Sj3] we detail the X-ray observations and the 
data analyses. Overall results and discussion are presented 
in <2] Throughout this paper, we adopt a cosmology with 
H Q = 70.5 km s- 1 Mpc~\ i2 M = 0.274, and Q. A = 0.726 (e.g., 
Komatsu et al. 2009). 

2. SAMPLE SELECTION 

We began our sample selection from the SDSS Data Re- 
lease 7 (DR7; Abazajian et al. 2009) quasar catalog (Schnei- 
der et al. 2010) which covers 9380 deg 2 of sky area. We 
searched for z > 4 RLQs satisfying one or both of the fol- 
lowing criteria: (1) log/? > 2.5 (R > 320); (2) M, < -29. 

9 Mi l defined "RLQs" as quasars with R > 100, while objects with 10 ^ 
R < 100 were classified as radio-intermediate quasars (RIQs). To simplify 
the terminology, we refer to all quasars with R ^ 10 as RLQs in this work. 
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Seven objects were first selected via this method. We also 
added another object (J0741+2520) which does not have 
SDSS spectroscopy, but was confirmed as a very luminous 
(Mi = -29.04) z = 5.194 RLQ in follow-up optical/UV spec- 
troscopy of radio-selected SDSS sources (McGreer et al. 
2009). Among the SDSS-covered z > 4 RLQs without sen- 
sitive X-ray coverage, these eight objects have the most re- 
markable radio and/or optical properties. They are ideal tar- 
gets, with high radio loudness, high optical luminosity, and 
high redshift, for economical Chandra snapshot observations. 
Six of the eight objects were awarded Chandra time in Cy- 
cle 12 (see Table [T] for an X-ray observation log). Four of 
the six targets (J1026+2542, J1412+0624, J1420+1205, and 
J1659+2101) have radio loudness log/? > 2.5. We define this 
radio-loudness value as the threshold for being "highly radio 
loud". RLQs satisfying this criterion represent the top ~ 5% 
of the RLQ population in terms of radio loudness (see the 
dash-dotted curve in Fig. [TJ. We note that Sbarrato et al. 
(2012b) have also recently studied J1026+2542, claiming 
that it is a blazar. The other two targets (J0741+2520 and 
J1639+4340) are only moderately radio loud (with log/? = 
1 .06 and log/? = 1 .62, respectively), but they are the most op- 
tically luminous (M, < —29) RLQs at z > 4 that have sensitive 
X-ray coverage (see Fig. [T}. 

In addition to our Chandra Cycle 12 targets, we further 
searched for other HRLQs at z > 4 satisfying our log/? crite- 
rion in the sky area north of 8 — —40° using the NASA/IPAC 
Extragalactic Database (NED)H The sky area north of 8 = 
—40° is covered by the 1.4 GHz NVSS survey, and thus we 
are able to calculate the radio loudness for each object. In 
fact, for typical z > 4 RLQs identified in current wide-field 
surveys (nij < 21), if an object satisfies our log/? > 2.5 cri- 
terion for being a HRLQ, it should have been detected by 
NVSS according to the NVSS sensitivity (~ 2.5 mJy). There- 
fore, our selection method should not introduce incomplete- 
ness owing to the radio-flux limitH For objects with both 
FIRST and NVSS detections, we adopt the radio-flux values 
from the FIRST catalog due to its better sensitivity; the ra- 
dio flux values in the FIRST and NVSS catalogs are generally 
consistent with each other. A total of 24 HRLQs at z > 4 were 
identified. We checked for sensitive archival X-ray coverage 
for these z > 4 HRLQs by Chandra, XMM-Newton, ROSAT\B 
or Swift. Thirteen have sensitive archival X-ray coverage (see 
Table QJEI while the other 11 do not (see TableEJlQ TableQ] 
and Table |2]together provide a complete list of known HRLQs 
at z > 4 in the sky area north of 8 = -40°. All of the 13 
objects having sensitive archival X-ray coverage are detected 

10 http://ned.ipac.caltech.edu/. 

1 1 However, it is still possible that some relevant HRLQs at z > 4 have not 
been identified due to the lack of optical/UV spectroscopic follow up. There- 
fore it is not possible to generate a fully complete sample of z > 4 HRLQs at 
present with available databases. 

12 We required the ROSAT observations to be pointed observations (i.e., 
not ROSAT All-Sky Survey) with an exposure time greater than 5 ks and an 
off-axis angle less than 19' (i.e., within the inner ring of the PSPC detector) to 
achieve adequate sensitivity. However, none of the z. > 4 HRLQs is covered 
by ROSAT observations satisfying these criteria. 

13 Only one known z > 4 HRLQ (PKS B 125 1-407) with sensitive X-ray 
coverage lies south of 5 = —40° , and thus is not included in our sample. See 
Yuan et al. (2006) for the X-ray properties of this object. 

14 The NED lists SDSS Jl 12429.62+283125.8 as a z. = 4.38 quasar based 
on SDSS DR6 data. However, the redshift of this object provided in the SDSS 
DR7 quasar catalog is z = 1.36. Therefore we do not consider this object as 
a z > 4 HRLQ although it satisfies our logR criterion. This object does not 
have sensitive X-ray coverage. 



in X-rays; all of them were targeted in their relevant X-ray 
observations except for SDSS J1235— 0003. The archival 
Chandra, XMM-Newton, or Swift observations of four objects 
(PMN Jl 155-3107, SDSS J1235-0003, GB 1713+2148, 
and PMN J 195 1+0 134) are reported in this work for the 
first time. Some of the archival objects have been classed 
as "blazars" and were covered by multiple X-ray observa- 
tions (see SJD. They were usually selected from radio (e.g., 
Green-Bank 6 cm, GB6, Gregory et al. 1996; the Parkes-MIT- 
NRAO survey, PMN, Griffith & Wright 1993) or X-ray (e.g., 
the ROSAT All Sky Survey, Voges et al. 1999) surveys. The 
individual X-ray properties of these blazars have been stud- 
ied in detail (see the last column of Table Q] for references). 
In this work, we combine the 13 archival sources and the 4 
newly observed HRLQs to form a sample to perform general 
systematic studies of HRLQs at high redshift. Our sample of 
HRLQs at z > 4 with sensitive X-ray coverage thus includes 
17 objects, eight of which have their X-ray properties pre- 
sented for the first time in this work. 

We searched for high-resolution radio images of these ob- 
jects obtained by the Very Long Baseline Array (VLBA) or 
Very Long Baseline Interferometry (VLBI); seven of them 
have available high-resolution images. Two objects have 
point-like source profiles and no indications of extended radio 
emission (J0913+5919, Momjian et al. 2004; J1235-0003, 
Momjian et al. 2004). Another three objects have clear ex- 
tended radio emission on milliarcsecond (mas) scales from 
relativistic jets although their extended radio emission only 
contributes < 10% of the total radio flux (Q 0906+6930, 
Romani et al. 2004; GB 1508+5714, Cheung 2004, Helm- 
boldt et al. 2007; GB 1428+4217, Cheung et al. 2012). 
The extended radio emission of GB 1508+5714 coincides 
with the extended X-ray emission, showing the signature of 
a relativistic jet (Cheung 2004). The extended radio feature 
of GB 1428+4217 in its VLA imaging at observed-frame 
1.4 GHz and 4.9 GHz (Cheung et al. 2012) lies 3.6" away 
from the core. It also coincides with the extended X-ray 
emission found in its Chandra observation. GB 1428+4217 
furthermore shows mas-scale extended radio emission in its 
VLBI images at observed-frame 2.3 GHz and 8.6 GHz. The 
extended radio emission of J1026+2542 is on a w 20 mas 
scale; it contributes w 40% of the total radio flux (Helmboldt 
et al. 2007). The mas-scale extended radio emission of the 
other object, J1659+2101, contributes w 30% of the total ra- 
dio flux at observed-frame 1 .6 GHz, while the image of this 
object at observed-frame 5 GHz shows no evidence of ex- 
tended radio emission (Frey et al. 2010)O The five objects 
with extended radio emission have higher radio loudness on 
average ((log/?) = 3.21 ±0.23) than the three objects without 
extended radio emission ((log/?) = 2.81 ±0.25). 

Besides the 17 HRLQs at z > 4, we will also separately dis- 
cuss the two highly optically luminous RLQs observed by 
Chandra in Cycle 12 even though they do not satisfy the 
"highly radio loud" criterion. All of our 19 objects are among 
the most-luminous RLQs in both the radio and optical/UV 
bands (see Fig. [2j. The monochromatic luminosity ranges of 

15 We also found high-resolution radio images for two z > 4 HRLQs 
without sensitive X-ray coverage listed in Table [2] The VLBI image of 
J08 13+3508 shows jet-like resolved structure at observed-frame 1.6 GHz, 
which contributes w 1 /4 of the total flux, while its image at observed-frame 
5 GHz does not have evidence of extended emission (Frey et al. 2010). 
J1242+5422 has extended radio emission in both of its VLBI images at 
observed-frame 1.6 GHz and 5 GHz, where the extended emission only con- 
tributes 3% of the total fluxes (Frey et al. 2010). 
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our objects are 10 33 -10 36 erg s 1 Hz 1 in the radio (at rest- 
frame 5 GHz) and 10 30 -10 33 erg s" 1 Hz" 1 in the optical/UV 
(at rest-frame 2500 A). GB 1713+2148 has remarkable prop- 
erties even compared to other HRLQs in our sample. It is 
much fainter (by a factor of s» 6) in the optical band, and it 
has the highest radio loudness (see Fig. [TJ. Its radio loudness 
is also higher than those of most of the objects (99%) in the 
full sample of Mil. All of our objects are spectroscopically 
confirmed broad-line quasars by either the SDSS or other ob- 
servations. They have precise redshift measurements based on 
their broad emission lines (e.g., Lya and/or C IV). The avail- 
able SDSS spectra of nine objects (five Chandra Cycle 12 
targets and four archival sources) are shown in Fig. |3] 

3. X-RAY DATA ANALYSIS 

The six RLQs at z > 4 targeted by Chandra in Cycle 12 
were observed with the S3 CCD of the Advanced CCD 
Imaging Spectrometer (ACIS; Garmire et al. 2003). The 
Chandra data were reduced using standard CIAO v4.3 rou- 
tines. We generated X-ray images for the observed-frame 
soft (0.5-2.0 keV), hard (2.0-8.0 keV), and full (0.5-8.0 keV) 
bands using ASCA grade 0, 2, 3, 4, and 6 events. The source 
detection was performed using the WAVDETECT algorithm 
(Freeman et al. 2002) with a detection threshold of 10~ 6 
and wavelet scales of 1, V2, 2, 2v / 2, and 4 pixels. All tar- 
gets were clearly detected by Chandra within 0.6" of the op- 
tical coordinates. To assess possible extended X-ray emission 
from a putative jet, we followed the method in §2.3 of Bassett 
et al. (2004) requiring a minimum of « 3 nearby counts off- 
set from the core by w 2" along roughly the same direction. 
None of our targets shows evidence for extended X-ray emis- 
sion. The X-ray counts were measured with aperture photom- 
etry using the IDL APER procedure. The aperture radius was 
adopted as 1.5" for each source (w 95% enclosed energy for 
the soft band, w 90% enclosed energy for the hard band; aper- 
ture corrections were applied). The background counts were 
retrieved from an annular region with inner and outer radii 
of twice and three times the aperture radius, and were scaled 
to the source-aperture area. All background regions are free 
of X-ray sources. Table |3]lists the X-ray counts in the three 
bands defined above, as well as the band ratio (defined as the 
ratio between hard-band and soft-band counts) and effective 
power-law photon index for each source. The effective pho- 
ton index was derived from the band ratio with the Chandra 
PIMMS3 tool, under the assumption of a power-law model 
with only Galactic absorption. 

The X-ray properties of the archival objects were obtained 
from the literature (see Table[T} except for those objects listed 
in Table [3] for which their archival Chandra, XMM-Newton, 
or Swift observations have not been previously published. 
For objects with multiple X-ray observations, we chose the 
ChandralXMM-Newton observation with the longest expo- 
sure time0 This approach should not introduce any biases 
in the following considerations regarding the X-ray emission 
strength. Vignali et al. (2003a) reported a tentative detection 
of GB 1713+2148 using aROSATHRl observation. This ob- 
ject was later observed with Chandra ACIS in Cycle 5 (PI: L. 
Gurvits). We processed the Chandra data for GB 1713+2148 
using the same method as described in the previous paragraph. 

1 6 http ://cxc . harvard.edu/toolkit/pimms .j sp 

17 We did not use the XMM-Newton observation with the longest exposure 
time for GB 1428+4217 because this object was undergoing an exceptional 
radio flare (see Worsley et al. 2006). 



The X-ray flux of GB 1713+2148 in the Chandra epoch is 
consistent with that in the ROSAT epoch, showing that this 
object did not have any strong X-ray variability. 

SDSS J1235— 0003 was serendipitously covered by an 
XMM-Newton observation on 2010 July 01. We performed 
data reduction and processing with standard XMM-Newton 
Science Analysis System (SAS; vlO. 0.0) routines. X-ray im- 
ages were generated for the observed-frame soft, hard, and 
full bands for the M0S1 and M0S2 detectors; this object is not 
covered by the pn detector. Source detection was carried out 
using the EBOXDETECT procedure. This object was detected 
only in the MDS1 full-band image at a w 3<T level, and thus 
we only used the M0 SI image for its X-ray photometry (it was 
not detected in the M0S2 images because it is very close to the 
CCD edge). We filtered the events file by removing the pe- 
riods of background flaring (49% of the total exposure time). 
The aperture radius for the photometry (53 .0") was the 90% 
enclosed-energy radius at 1 .5 keV based on the point spread 
function (PSF) of the M0S1 detector at an off-axis angle of 
12.8'. The upper limits on X-ray counts in the soft and hard 
bands were calculated as 3^/N, where N is the total counts in 
the aperture. The vignetting was corrected by calculating the 
effective exposure time (4.2 ks) at the source location based 
on the exposure map. 

Another two archival objects, PMN Jl 155— 3107 and 
PMN J 195 1+0 134, were targeted by the Swift X-ray Tele- 
scope (XRT; Burrows et al. 2005) in two (ObsIDs: 
00036263001/2, total exposure time 5.1 ks) and three (Ob- 
sIDs: 00036791001/2/3, total exposure time 10.4 ks) obser- 
vations, respectively. For each object, we first generated the 
exposure maps for all XRT events files using the XRTEX- 
POMAP routine. We then merged the event files and expo- 
sure maps for each object using the XSELECT and XIMAGE 
packages, respectively. Similar to the approach for the Chan- 
dra and XMM-Newton data, XRT images were created for the 
observed-frame soft, hard, and full bands. Source detection 
was performed with the WAVDETECT algorithm with a de- 
tection threshold of 10~ 6 . PMN Jl 155-3107 is detected in 
the full band and soft band, while PMN J 195 1+0 134 is de- 
tected in all three bands. Aperture photometry was performed 
using aperture radii of 59.6" and 64.9", which are the 90% 
enclosed-energy radii in the full band based on the PSF of 
the XRT detector at off-axis ang les of 1.0' and2.6'0 respec- 
tively. The X-ray photometry of these two archival objects is 
also listed in Table [3] 

Table |4]lists the key X-ray, optical, and radio properties of 
our sample: 

Column (1): The name of the quasar. 

Column (2): The apparent SDSS /-band magnitude of the 
quasar. For objects covered by the SDSS footprint, the val- 
ues are obtained either from the SDSS DR7 quasar catalog 
(BEST photometry; Schneider et al. 2010) or from the SDSS 
database. For other objects, the ;-band magnitude was con- 
verted from AB 1450(1+2) under the assumption of a power-law 
spectral index of a v — —0.5 (/ v « V av ; e.g., Vanden Berk 
et al. 2001). The values of ABi45Q(i +z ) were calculated from 
/?-band magnitudes (from the NED) using the empirical rela- 
tion of AB 145 o(i +z ) = R - 0.684z+ 3.10 (Kaspi et al. 2000; 
Vignali et al. 2003a). 

18 The Swift mission has excellent rapid pointing capability, while its point- 
ing accuracy is not exceptional. Therefore, it is normal to have off-axis angles 
of sa 1'— 2' for targeted observations. 



HIGHLY RADIO-LOUD QUASARS ATz > 4 



5 



Column (3): The absolute SDSS /-band magnitude for the 
quasar, M,, from the SDSS DR7 quasar catalog, which was 
calculated from the apparent SDSS /-band magnitude in Col- 
umn 2 by correcting for Galactic extinction and assuming a 
power-law spectral index of a v = —0.5. 
Column (4): The Galactic neutral hydrogen column den- 
sity calculated with the Chandra C0LDEN3 tool (Dickey & 
Lockman 1990; Stark et al. 1992), in units of 10 20 cm" 2 . 
Column (5): The count rate in the observed-frame soft X-ray 
band (0.5-2.0 keV) for the Chandra observed objects, in units 
of 10~ 3 s _1 . 

Column (6): The X-ray flux in the observed-frame soft band 
(0.5-2.0 keV) corrected for Galactic absorption. This mea- 
surement was obtained with the Chandra PIMMS tool and 
is in units of 10~ 14 erg cm~ 2 s _1 . An absorbed power-law 
model was used with the photon index (Fx) listed in Column 9 
and the Galactic neutral hydrogen column density (Nh) listed 
in Column 4. For objects without Tx information or only with 
lower limits upon Tx, we adopt Tx = 1 .6 which is typical for 
radio-loud quasars (e.g., Page et al. 2005). 
Column (7): The X-ray flux density at rest-frame 2 keV 
generally obtained from the count rate in the observed- 
frame 0.5-2.0 keV band with PIMMS, in units of 10~ 32 
erg cm~ 2 s _1 Hz -1 and corrected for Galactic absorption. 
Although for our objects rest-frame 2 keV corresponds to 
observed-frame 0.3-0.4 keV, we chose the standard approach 
to using the observed-frame 0.5-2.0 keV band because we are 
able to minimize effects from potential intrinsic X-ray absorp- 
tion and X-ray instrumental contamination below observed- 
frame 0.5 keV, and we can also obtain better photon statistics 
with the larger numbers of X-ray counts in this band. 
Column (8): The logarithm of the quasar luminosity (erg s ) 
in the rest-frame 2-10 keV band corrected for Galactic ab- 
sorption. 

Column (9): The power-law photon index (Fx) of the X- 
ray spectrum. For objects listed in Table [3] the values were 
obtained from band-ratio analysis (see the last column of 
Table [3}. For the other objects, the values were retrieved 
from the literature. Some objects have shown X-ray spec- 
tral variability (e.g., GB 1428+4217; see Fig. 4 of Wors- 
ley et al. 2006). The choice of Tx values do not signifi- 
cantly affect the X-ray flux calculation. We have chosen the 
Chandra/XMM-Newton observations with the longest expo- 
sure time which should best constrain the power-law photon 
index. 

Column (10): The continuum flux density at rest-frame 
2500 A in units of 10~ 27 erg cm~ 2 s _I Hz -1 . For objects 
appearing in the SDSS DR7 quasar catalog, the values were 
obtained from the catalog of Shen et al. (2011). For other 
SDSS covered objects, the flux density was calculated using 
the SDSS photometry with the composite spectrum for SDSS 
quasars in Vanden Berk et al. (2001). For objects not covered 
by the SDSS, the flux density was calculated using the opti- 
cal photometry from the literature under the assumption of a 
power-law spectral index of a v = —0.5. 
Column (11): The logarithm of the monochromatic luminos- 
ity (erg s _I Hz -1 ) at rest-frame 2500 A, which was calcu- 
lated from the flux density at rest-frame 2500 A given in Col- 
umn 9. A standard cosmological bandpass correction was ap- 
plied when converting flux into monochromatic luminosity. 
Column (12): The radio power-law slope Or (f v V a ') be- 

19 http://cxc.harvard.edu/toolkit/colden.jsp 



tween observed-frame 1.4 GHz and 5 GHz. The 1.4 GHz 
flux density was obtained from the FIRST or NVSS surveys, 
while the 5 GHz flux density was obtained from the GB6 sur- 
vey, the PMN survey, or from the literature for individual ob - 
jects (J1639+4340, Holt et al. 2004; Q 0906+6930, Romani 
2006; J0913+5919, Momjian et al. 2004). The observations 
at 1.4 GHz and at 5 GHz are not simultaneous. All objects 
with Or measurements are flat-spectrum radio quasars (FS- 
RQs; Or > -0.5) except for J0913+5919. For the five objects 
without 5 GHz flux measurements, we assume = to keep 
consistency with other objects in our sample. 
Column (13): The logarithm of the monochromatic luminos- 
ity (erg s ~ 1 Hz 1 ) at rest-frame 5 GHz, obtained from the flux 
density at rest-frame 5 GHz, f$ ghz- This quantity was calcu- 
lated using a radio power-law slope given in Column 12 and 
the flux density at observed-frame 1.4 GHz, /1.4 ghz, which 
were obtained from the FIRST or the NVSS surveys. A stan- 
dard cosmological bandpass correction was also applied. 
Column (14): The logarithm of the radio-loudness parameter, 
given by 

io g *=iog(|^y a) 

V / 4400 A / 

The denominator, /4400 a> was calculated from ^9500 A us ~ 
ing an optical power-law slope of a v = —0.5. The numera- 
tor, fs GHz, was calculated using the same method as in Col- 
umn 13. All measures of flux density are per unit frequency. 
The definition of radio loudness in Mil used /2500 A as the 
denominator. We converted their values to our definition for 
the scientific analyses described in the following section. 
Column (15): The a ox parameter, defined by 



anx= 1Og(/2kev//25Q0A) =0-384 log ( /2keV 



log(V 2 kev/v 25 ooA) 



2500 A 



(2) 



which represents the slope of an assumed power law connect- 
ing rest-frame 2500 A and 2 keV. Our UV and X-ray mea- 
surements were not simultaneous, and thus the a ox parameter 
could be affected by variability. 
Column (16): Aa ox .RQQ, defined as 

AOo^rqq = «ox (measured) - a ox ,RQQ (expected). (3) 

The expected CV.rqq value for a typical RQQ is calculated 
from the ao X -L 2 5oo A correlation given as Equation (3) of Just 
et al. (2007). The statistical significance of this difference 
compared to RQQs (given in parentheses) is in units of <j, 
which is given in Table 5 of Steffen et al. (2006) as the RMS 
of otox for several ranges of luminosity. The Aa ox .RQQ pa- 
rameter quantifies the excess of X-ray emission from the rel- 
ativistic jet of the RLQs compared to that of RQQs for which 
the X-ray e miss ion is mainly from the accretion disk and its 
corona. In ! j4.2l we will validate the use of Aa 0X ,RQQ as a use- 
ful diagnostic quantity for HRLQs. 
Column (17): Aa ox .RLQ, defined as 

A0£ OX! r L q = a ox (measured) - a 0Xj R L Q (expected) . (4) 

The expected a ox ,RLQ value for a typical RLQ is obtained 
from the L 2 keV-^ 2 5oo A~^5 GHz correlation based on the full 
sample in Ml 1 given in their Table 7, which is equivalent to 

cv,rlq = -0.199 L 2500 A + 0.105 L 5 GHz + 1.194. (5) 

The Aa ox .RLQ parameter assesses the relative X-ray brightness 
of a RLQ compared to typical RLQs (mostly at z = 0.3-2.5 
and with a median redshift of z — 1 .4) with similar optical/UV 
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and radio luminosities that have both disk/corona-linked and 
jet-linked X-ray emission. 

4. RESULTS AND DISCUSSION 
4. 1 . Enhanced X-ray Emission at High Redshift 
4.1.1. The Relative X-ray Brightness of z > 4 HRLQs 

To assess the enhancement of the X-ray emission of our 
objects at z > 4 compared to similar objects at lower red- 
shift, we compare their distributions of a ox , Aa ox rqq, and 
AOox.rlq to those of RLQs at z < 4. Fig. [4] shows the po- 
sitions of our objects and the typical RLQs in Mil in the 
logfl-(a ox , AOo Xi rqq, Ao£ox,rlq) planes. Our objects gen- 
erally occupy the same region as the typical RLQs in the 
logR-a ox plane (see the top panel of Fig. [4]). However, our 
sample appears to have an excess of objects with stronger 
X-ray emission in the logR-(Aa ox ,rqq, AOox.rlq) planes. 
The thick black lines in the middle panel of Fig. |4]show the 
mean Ac<o X ,rqq values of Ml 1 objects binned by radio loud- 
ness (Alog/? = 0.2 per bin). Fifteen out of our 17 HRLQs 
at z > 4 have Aa ox rqq values greater than the mean value of 
the corresponding log/? bin. Similarly, 14 of our 17 HRLQs 
at z > 4 have AOo X ,rlq values greater than zero, which gener- 
ally represents the mean Aoc ox .rlq value of Ml 1 objects (see 
the bottom panel of Fig. |4j». We will show later that the gener- 
ally higher Aa ox .RQQ and Aa 0X ,RLQ values of our objects are 
caused by stronger X-ray emission rather than weaker opti- 
cal/UV emis sion (see the broad-band SED studies on our ob- 
jects in ^4.41 . It is worth noting that the X-ray weakest object 
in our z > 4 HRLQ sample (J0913+5919; see the red square 
with the smallest Aa ox rqq and Aa ox rlq values in Fig.[4]i has 
a steep radio spectral slope (0^- = —0.67), while the others 
with available c^- values are all FSRQs. 

To verify this excess quantitatively, we obtained the mean 
values and statistical distributions of Ao! ox .rqq and Aa ox .RLQ 
for our HRLQs and a comparison sample consisting of ob- 
jects in the full sampl^] of Ml 1 with \ogR > 2.5 and z < 4. 
Our sample, which includes objects originally selected in the 
optical (7 objects; e.g., SDSS), radio (9 objects; e.g., GB6 
and PMN), or X-ray (1 object; from the ROSAT All-Sky Sur- 
vey) bands, is somewhat heterogeneous. To minimize poten- 
tial selection biases (e.g., the most exceptional objects be- 
ing preferentially targeted in X-ray observations), we need a 
sample containing most of the known z > 4 HRLQs regard- 
less of their selection methods and having nearly complete 
X-ray coverage. To achieve this goal, we chose an optical- 
magnitude cut of m, < 20 on all known z > 4 HRLQs north 
of 8 = -40° (see Tableland TableEJ. There are 15 HRLQs 
at z > 4 satisfying this m, cut, and 12 of them have sensitive 
X-ray coverage (and thus are included in our sample). For 
the three HRLQs without sensitive X-ray coverage, two of 
them (J0813+3508 and J1242+5422, both optically selected) 
were in our original Chandra Cycle 12 proposal, but were not 
awarded Chandra time (see jj2). These two objects have the 
highest redshifts, optical brightnesses, and optical luminosi- 
ties among those z > 4 HRLQs without sensitive X-ray cov- 
erage (see Tabled. The /-band magnitude of the other HRLQ 
(J2220+0025, optically selected) is close to our m, cut. To 
maintain consistency, we also applied the m, < 20 criterion to 
our comparison sample from Mil. There are 283 objects in 
total in our comparison sample, including 269 X-ray detected 

20 The full sample of Ml 1 contains a mixture of objects that are optically 
selected (~ 80%), radio selected (~ 15%), or X-ray selected (~ 5%). 



objects and 14 objects with X-ray upper limits. 

The mean values of Aa 0X ,RQQ and Aa ox .RLQ shown in Ta- 
ble [5] were calculated using the Kaplan-Meier estimatoiFI 
implemented in the Astronomy Survival Analysis (ASURV) 
package (e.g., Lavalley et al. 1992). These mean values indi- 
cate that HRLQs at z > 4 have stronger X-ray emission (by a 
factor of w 3 on average) than those at z < 4 with similar UV 
and/or radio luminosity. The distributions of Aa ox rqq and 
Acv rlq are shown in Fig. [5] The Peto-Prentice test (e.g., 
Latta 1981), also implemented in the ASURV package, was 
employed to determine whether the z > 4 and z < 4 HRLQs 
follow the same Aa ox .RQQ and AOo X ,rlq distributions. The 
Peto-Prentice test is preferred over other two-sample tests be- 
cause it is the least affected by factors such as unequal sample 
sizes (e.g., Latta 1981) which do exist in our study. The dis- 
tribution of Aa ox .RQQ for our objects is significantly different 
from that of the HRLQs at z < 4. The null-hypothesis (i.e., 
the two samples following the same distribution) probability 
is only 1.21 x 10~ 4 . The two-sample test for the Aa ox .RLQ 
distribution reaches a similar conclusion; the null-hypothesis 
probability is slightly higher (1.37 x 10~ 4 ), but it still shows 
that the two samples follow different distributions at a > 3(7 
level. One might speculate that the significant differences of 
the Aa ox distributions are caused by a "tail" of z > 4 objects 
with extraordinary X-ray brightness. However, Fig. |5] sug- 
gests that the entire AOo X histograms of our z > 4 HRLQs 
are shifted toward high Aa ox values relative to those of z < 4 
HRLQs (also see Fig. |4j. Therefore, our sample of HRLQs 
shows a significant X-ray emission enhancement at z > 4 
over those HRLQs at lower redshift. However, one should 
keep in mind that although we have compiled a sample of 
Z > 4 HRLQs with nearly complete X-ray coverage using all 
available databases, potential selection biases could still exist. 
For example, although we have assembled all known z > 4 
HRLQs north of 8 = —40°, this is still not a fully complete 
sample (see Footnote 11). Therefore, in the next subsection 
we describe a variety of tests of the robustness of our results. 

4.1.2. Robustness Tests 

In order to test the robustness of our results against selection 
issues, we carried out further two-sample analyses by con- 
structing other comparison samples or controlling for a vari- 
ety of relevant parameters (see results in Table |5): 

1 . Many of our z > 4 HRLQs are radio selected, while the 
majority of our comparison sample from Mil is op- 
tically selected. To minimize the effects of different 
selection methods, we first compare our radio-selected 
z > 4 HRLQs with a radio-selected subsample in Mil. 
We chose their "Einstein" supplemental sample (see 
§2.2.1 ofMll) which was primarily radio selected. The 
X-ray observations of these objects were first presented 
in Worrall et al. (1987), and we transform their lumi- 
nosities to our adopted cosmology. We also incorporate 
improved X-ray coverage for a handful of objects, as 
in Ml 1. All of these objects have z < 4 and m, < 200 
we further require them to be HRLQs (i.e., \ogR > 2.5). 
Our z > 4 HRLQs have significantly different Aa ox dis- 
tributions (both Aa ox .RQQ and Aoc ox rlq) from those of 

21 The Kaplan-Meier estimator is applicable to censored data; note some 
of the Mil data points are censored. 

22 The m, values were estimated from their monochromatic luminosities at 
rest-frame 2500 A. 
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this radio-selected comparison sample at a ps 3d level. 

2. We generate a quasi-radio-selected comparison sample 
following §2.1.1 of Mil by only selecting those ob- 
jects that were targeted by SDSS spectroscopy owing 
to being FIRST radio sources (i.e., with the "FIRST" 
flag set in the SDSS quasar catalog)0 Our z > 4 
HRLQs have significantly different Aa ox distributions 
from those of this "radio-selected" comparison sample 
at a ps 4a level. 

3. As discussed in S|2] our z > 4 HRLQs are among the 
RLQs with the highest radio luminosities; this may in- 
troduce a luminosity bias into our sample. As discussed 
in Ghisellini et al. (2011), for example, more lumi- 
nous blazars have relatively more hard X-ray emission 
(because the peak photon energy shifts to lower fre- 
quencies), and this could make the observed X-ray flux 
greater for high-redshift blazars. Therefore, we con- 
struct another comparison sample from Mil to have 
comparable radio luminosity (logLs ghz > 34; compare 
with Fig. |2). After controlling for radio luminosity, 
the significance of the difference in the Aa ox distribu- 
tions becomes slightly lower, but still remains at a ps 3(7 
level. 

4. Similarly, we construct another comparison sample to 
have comparable optical/UV luminosity (log L2500 A > 
30.9; compare with Fig.|2]i to our z > 4 HRLQs to avoid 
optical/UV luminosity bias. The significance of the 
Aa ox distribution difference is at a > 2.5(7 level after 
controlling for optical/UV luminosity. 

5. All of our z > 4 HRLQs with available Or are FS- 
RQs {Ok > -0.5; see © except for JO913+591S0, 
which may suggest that low-inclination RLQs are over- 
represented in our sample. Therefore, for another com- 
parison sample, we selected the Mil full-sample ob- 
jects with available measurements that satisfy the 
FSRQ criterion. After controlling for this effect, our 
results still have > 3.5(7 significance. 

In summary, our finding that z > 4 HRLQs have a substantial 
X-ray emission enhancement over HRLQs at lower redshifts 
remains significant after controlling for those factors which 
may introduce biases. We also note that, owing to the re- 
spectable sizes of the samples being compared, variability ef- 
fects upon Aa ox distributions should tend to average out (also 
see §3 for explanation of how our selection of the X-ray ob- 
servations utilized was designed to avoid biases due to vari- 
ability). 

It is possible that this X-ray enhancement of HRLQs is not 
solely confined to z > 4 but arises at lower redshift. To inves- 
tigate this briefly, we compared the Aa ox rqq and Ao! ox .rlq 
distributions of z = 3-4 HRLQs (9 objects) in the full sam- 
ple of Mil with those of z < 3 HRLQs (274 objects), and 
found they are different at a >5(7 level (see Table [6j z = 3-4 
HRLQs have both larger Aa ox .RQQ and Aa ox .RLQ values which 
indicate an X-ray emission enhancement by a factor of ps 3). 
However, the Aoc ox .rqq values of z = 2-3 HRLQs (39 objects) 

23 This is not a true radio-selected sample because, e.g., they were re- 
stricted by SDSS magnitude limits for FIRST sources. See §2. 1 . 1 of M 1 1 for 
details. 

24 This object is not included in our two-sample tests since it does not 
satisfy m, < 20. 



follow a similar distribution to that of z < 2 HRLQs (235 
objects; 0.18(7 level difference). For the Ao! ox ,rlq parame- 
ter, although its distribution for z — 2-3 HRLQs is different 
from that of z < 2 HRLQs at a w 2(7 level, the mean Aa ox .RLQ 
value of z = 2-3 HRLQs is slightly lower than that of z < 2 
HRLQs (-0.030 ±0.019 vs. 0.005 ± 0.008). Therefore, the 
above tests suggest that the X-ray emission enhancement of 
HRLQs begins to arise atz ps 3. This statistically independent 
result from z = 3-4 also supports the robustness of our find- 
ings of a HRLQ X-ray emission enhancement at z > 4. Fur- 
thermore, after combining the samples of z > 4 HRLQs and 
z = 3-4 HRLQs, we find that z > 3 HRLQs have an X-ray 
emission enhancement over z < 3 HRLQs at a > 5.7(7 signif- 
icance level corresponding to a null-hypothesis probability of 
ps 10~ 8 (also see Table©. 

4. 1 .3. Relevant Physical Considerations Regarding the X-ray 
Emission Enhancement 

Our results suggest that the most-luminous relativistic jets 
generate stronger X-ray emission in the early universe. Previ- 
ous studies by Bassett et al. (2004) and Lopez et al. (2006) that 
included a mix of HRLQs (including a subset of those in our 
sample) and moderately radio-loud quasars found that high- 
redshift objects did not show substantially stronger X-ray 
emission compared to a lower-redshift sample of RLQs de- 
tected by ROSAT (from Brinkmann et al. 1997). We as- 
sess further these previous results by performing a similar 
two-sample analysis on the moderately radio-loud quasars 
(1 < log/? < 2.5) in the full sample of Mil (which includes 
the moderately radio-loud objects from Bassett et al. 2004 and 
Lopez et al. 2006). The Acv.rqq distribution of z > 4 moder- 
ately radio-loud quasars (10 objects; all X-ray detected) is dif- 
ferent from that of z < 4 moderately radio-loud quasars (280 
objects with 48 X-ray upper limits) at a ps 2(7 level, while 
the Aa ox rlq values of z > 4 moderately radio-loud quasars 
follow a similar distribution to that of z < 4 moderately radio- 
loud quasars (0.57(7 level difference). Therefore, there is no 
firm evidence that moderately radio-loud quasars have a sig- 
nificant X-ray emission enhancement compared to those at 
lower redshift (recall also from §Q] that RQQs show no evo- 
lution in X-ray properties with redshift). This suggests that 
the X-ray enhancement in high-redshift HRLQs occurs within 
the jet-linked component, which likely (e.g., Mil) contributes 
an increasing fraction of the X-ray continuum with increasing 
radio-loudness value sF^I 

The enhancements in relative nuclear X-ray emission we 
find for z > 4 HRLQs are less than would be predicted for 
an X-ray jet-linked component dominated by IC/CMB emis- 
sion, which possesses a strong (1 +z) 4 dependence on red- 
shift (e.g., Schwartz 2002; see §1). For example, the median 
redshift for the RLQs from Miller et al. (2011) with z < 4, 
log/? > 2.5, and m.\ < 20 is z = 1.3. If the jet-linked contribu- 
tion to the nuclear X-ray emission is at least ps 50% at z = 1.3, 
as is a conservative estimate for RLQs with log/? > 2.5 (e.g., 
Zamorani et al. 1981; Figure 7 of Mil), then for a pure 
IC/CMB X-ray jet the overall enhancement in the nuclear 
X-ray emission at z = 4.4 would be >16, corresponding to 
an increase in Aa ox rlq of 0.46. However, our z > 4 HRLQs 
only show an increase in A oc ox rlq of 0.15-0.20 (see Ta- 
ble 5). Instead, our results are consistent with an X-ray jet- 

25 It is obviously also possible, but not required, that the jet-linked X-ray 
generation mechanism of HRLQs differs from (and is more sensitive to red- 
shift than) that of moderately radio-loud quasars. 
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linked component in which the fractional contribution from 
the IC/CMB process at z = 1.3 is <6% of the nuclear X-ray 
emission, with the balance of the jet-linked X-ray emission 
arising from inverse Compton scattering of other seed pho- 
ton fields such as radiation from the broad-line region or dust 
(e.g., Sikora et al. 2009), from the jet synchrotron emission 
itself (e.g., Sokolov et al. 2004; Meyer et al. 2012 find that 
synchrotron self-Compton emission is relatively less impor- 
tant in high-powered jets), or from non-cospatial synchrotron 
radiation in a structured jet (e.g., Migliori et al. 2012). The 
overall enhancement to the nuclear X-ray continuum then 
rises more gradually with redshift; at z = 4.4/3.0/2.0 RLQs 
would be X-ray brighter by factors of <2. 8/1. 5/1.1 (relative 
to a z = 1.3 template), corresponding to Aa ox rlq increases 
ofO.17/0.07/0.02. 

Another possible explanation for the X-ray enhancement is 
that an increasing contribution of the X-ray emission from the 
jets is due to the inverse-Compton scattering of the photon 
field of the host galaxy; e.g., Hardcastle & Croston (2011) 
provided a detailed modeling of upscattering of host-galaxy 
photon fields into high-energy emission. The X-ray en- 
hancement of our sample can be explained if the host galax- 
ies of high-redshift HRLQs are more luminous in the in- 
frared/optical band than those at lower redshift. This higher 
infrared/optical luminosity could be generated by the en- 
hanced star-formation activity of high-redshift host galaxies 
(e.g., Wang et al. 2008, 2011a; Mor et al. 2012). Smail 
et al. (2012) found that for the high-redshift (z =3.6-3.8) 
radio galaxies in their study, the luminosity of the extended 
X-ray halo is correlated with the far-infrared luminosity of 
the galaxy, which supports the X-ray generation mechanism 
of inverse Compton scattering of the far-infrared photons pro- 
duced by star-formation activity. Since the cosmic star for- 
mation rate density has a broad plateau starting at z ~ 1 and 
running to z ~ 4 (e.g., Hopkins & Beacom 2006), one might 
expect in this model the enhancement to emerge at z ~ 1 rather 
than at z ^3-4. However, the generally massive host galax- 
ies of quasars may show stronger evolution effects at earlier 
cosmic times owing to cosmic downsizing (e.g. Cowie et al. 
1996; Panter et al. 2007). For example, Archibald et al. 
(2001) found a (1 +z) 3 evolution of the total infrared lumi- 
nosity of radio galaxies at high redshifts. 

The a ox and AOo X parameters are defined based on the 
X-ray flux density at rest-frame 2 keV, which corresponds to 
observed-frame ~ 0.3-0.4 keV for our z > 4 HRLQs. How- 
ever, the rest-frame 2 keV flux density was derived from 
the observed-frame 0.5-2.0 keV count rate (see Sj3]l, which 
means that we have extrapolated the X-ray spectrum to rest- 
frame 2 keV assuming a single power-law model. It is known 
that some z > 4 HRLQs show X-ray spectral flattening below 
observed-frame 0.5 keV, i.e., an X-ray flux deficit compared 
to a single power-law model (e.g., Yuan et al. 2006). There- 
fore, the extrapolation process may have over-estimated the 
X-ray flux density at rest-frame 2 keV. Although X-ray spec- 
tral flattening has also been found in a few z < 3 objects (e.g., 
Page et al. 2005), it would not affect the X-ray flux calcu- 
lation because rest-frame 2 keV is covered by the observed- 
frame 0.5-2.0 keV band for these objects. To test whether the 
X-ray emission enhancement we found for z > 4 HRLQs is 
caused by this kind of X-ray spectral curvature, we recalcu- 
lated the a ox and Aa ox values for our z > 4 HRLQs utilized 
in the two-sample analyses above (i.e., those with m, < 20) 
from the count rates in the observed-frame ultrasoft band 



(0.3-1.0 keV; see results in Table |7]i except for the two Swift 
sources (PMN Jl 155-3107 and PMN J1951+0134) which 
do not have adequate counts in the ultrasoft band due to the 
limited sensitivity of the Swift XRT These a x,us values ob- 
tained from the ultrasoft band are similar to those we have ob- 
tained in |j3j the mean difference between them is (a OXj us — 
a ox ) = —0.03 ±0.01. If we use the Aa ox values calculated 
from the ultrasoft band (Aa ox rqq us an d Aa ox .rlq.us) m 
the two-sample tests above, the AOo X distributions of z > 4 
HRLQs and z < 4 HRLQs are still different at a > 3(7 level. 
Therefore, X-ray spectral curvature only contributes to a small 
fraction, if any, of the X-ray emission enhancement of z > 4 
HRLQs over those at lower redshift. In fact, some of the z > 4 
HRLQs without the apparent X-ray spectral curvature (e.g., 
GB 1508+5714; Yuan et al. 2006) are among the objects we 
find to have the strongest X-ray enhancements. 

Considering the seven z > 4 HRLQs with high-resolution 
radio images (see Sj2), the five with extended radio 
emission on mas scales have higher AOox.rqq values 
((AcV rqq) = 0.51 ±0.09) on average than those of the two 
without extended radio emission ((Aa ox rqq) = 0.15 ± 0.24). 
This might be because the objects with extended radio emis- 
sion have higher average radio loudness (see $2). How- 
ever, After considering the strength of their radio emis- 
sion, the objects with extended radio emission still have a 
higher average Ag^rlq value ((Aa ox ,RLQ) = 0.23 ± 0.09) 
than that of the objects without extended radio emission 
((Aa ox , rlq) — —0.07 ±0.23) although with substantial un- 
certainty due to the limited number of objects. 

4.2. Optical/UV Broad Emission Lines 

The optical/UV spectra of quasars at z > 4 usually show 
prominent broad emission lines such as Lya + N V and C IV 
(see Fig. [3). For RLQs, if the optical/UV continuum emission 
is also beamed following the radio emission, one might expect 
that the emission-line rest-frame equivalent widths (REWs) 
would be correlated with the beaming of the radio emission. 
Therefore, we investigate the relation between emission-line 
REWs and radio loudness for our sample of HRLQs, and 
compare this relation to that for the majority of RLQs over 
a wide range of redshift. 

For our z > 4 HRLQs with SDSS spectroscopy, we ob- 
tained measurements of REWs from Diamond-Stanic et al. 
(2009) and Shen et al. (2011) for Lya + N V and C IV, 
respectively (see Table HJ0 Two objects (J1026+2542 
and J1420±1205) with SDSS spectroscopy do not have 
REW(Lya + N V) measurements in Diamond-Stanic et al. 
(2009)0 We measured REW(Lya + N V) for these two 
sources following the method in §2 of Diamond-Stanic et al. 
(2009). Available REW measurements for other objects with- 
out SDSS spectroscopy were also obtained from the literature 
(see Table [8]). Measurements of REW(Lya + N v) are usu- 
ally affected by the presence of the Lya forest, so our values 
likely underestimate the intrinsic strength of Lya. Although 
the number density of Lya absorbers is known to evolve with 
redshift (e.g., Weymann et al. 1998; Songaila & Cowie 2002; 
Janknecht et al. 2006), this effect will not introduce biases 
to our correlation analyses of emission-line REWs with radio 

26 We also list the values of J1639+4340 in Tableland the following 
Table|9]for reference despite it is a moderately radio-loud quasar. 

27 Diamond-Stanic et al. (2009) measured the REW(Lya + N V) values 
for SDSS DR5 quasars (Schneider et al. 2007) at z > 3. These two objects 
are in the SDSS DR7 quasar catalog but not in the DR5 quasar catalog. 
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loudness assuming the RLQ radio-loudness distribution does 
not evolve with redshift. The number of objects in our sample 
that have emission-line REW measurements is not adequate 
for reliable correlation analysis. Therefore, we first investi- 
gate these correlations for a sample of RLQs across a wide 
range of redshift, and then see whether our z > 4 HRLQs fol- 
low similar trends. We selected all SDSS DR7 RLQs which 
have REW(Lya + N V) and REW(C IV) measurements from 
Diamond-Stanic et al. (2009) and Shen et al. (2011), respec- 
tively. Fig.[6]shows the relation between emission-line REWs 
and radio loudness, for which highly significant positive cor- 
relations are found (see Table [TOlfor the correlation probabil- 
ities). Fig. [6] also shows the best-fit linear correlation in the 
logarithmic parameter space using the IDL LINFIT procedure, 

loaREW(Lya + N V) = (1 .68 ± 0.06) + (0.05 ± 0.02) logfl, 

(6) 

and 

logREW(C IV) = (1.20±0.02) + (0.16±0.01)logfl. (7) 

Our HRLQs generally follow the above positive correla- 
tions between emission-line REWs and radio-loudness (see 
the filled stars and squares in Fig. [6), showing that these cor- 
relations do not seem to evolve with redshift. These pos- 
itive correlations are consistent with the results in Kimball 
et al. (201 la) on REW(Mg II) and REW(C IV) for SDSS DR5 
quasars0 One possible explanation for them is that the broad 
emission-line region of RLQs may emit anisotropically (see 
§6.2 of Kimball et al. 2011a). Additionally or alternatively, 
the decrease in synchrotron peak frequency with increasing 
peak or kinetic luminosity in FSRQs (e.g., Meyer et al. 201 1, 
their Figure 4) might more stringently limit jet dilution of op- 
tical emission-line features in more radio-loud quasars. Our 
results indicate that for HRLQs with broad emission lines 
(i.e., not BL Lac objects), the beamed relativistic jets gen- 
erally have negligible contribution to the optical/UV continua 
(except for a few outliers with exceptionally weak optical/UV 
emission lines; see below). This is consistent with recent re- 
sults from SED studies of HRLQs (e.g., Ghisellini et al. 2010) 
that their optical/UV emission is disk dominated. Further- 
more, this finding of a minimal jet-linked contribution to the 
optical/UV continuum validates the use of Acv.rqq as a di- 
agnostic quantity for our HRLQs; note in Column (16) of Ta- 
ble |4]that the a O x-^2500 A correlation for RQQs is utilized. 

To assess further the jet contribution to the optical/UV con- 
tinua of our HRLQs, we calculated the non-thermal domi- 
nance parameter (NTD = L bs/ipred) of the optical/UV con- 
tinuum following Shaw et al. (2012). L Q b s is the observed 
monochromatic continuum luminosity near an emission line 
(VL1350 near C IV in our case), while Lpnxi is the predicted 
VL1350 from the relation between continuum luminosity and 
emission-line luminosity in Shen et al. (2011; see their 
Equation 14) obtained from typical quasars. Lp re d represents 
the optical/UV continuum emission from the accretion disk. 
NTD^> 1 means that there is substantial jet-contributed op- 
tical/UV continuum. We calculated the NTD values for our 

28 Kimball et al. (2011a) found positive correlations between the core 
radio-to-optical flux ratio Ri and the REW values of Mg II and C IV for SDSS 
RLQs. The definition of Ri is generally different from that of radio loudness 
because the latter is defined based on the total radio luminosity. However, 
the positive correlations in Kimball et al. (2011a) also hold true for their 
subsample of RLQs without extended radio emission (the "core" morphology 
class which makes up s« 70% of their full sample; see their Table 8). In this 
case, the definition of Ri is similar to that of radio loudness. Therefore, our 
results are consistent with those of Kimball et al. (201 la). 



HRLQs with SDSS spectroscopy covering the C IV line (see 
Table |9]). Three objects (J1420+1205, J1235-0003, and 
Jl 325+1 123) have NTD values close to unity, supporting our 
results that the relativistic jets of our HRLQs generally have 
little contribution to their optical/UV continua. Another ob- 
ject, GB 1508+5714 (NTD = 1.3), may possibly have some 
jet-linked optical/UV continuum emission. However, its NTD 
value is consistent with unity considering the scatter of the 
relation between VL1350 and Lq rv (~ 0.2 dex; see §3.7 of 
Shen et al. 2011). The other two HRLQs (J1412+0624 and 
J1659+2101) have NTDs ss 6 due to their exceptionally weak 
C IV lines that are discussed further below. 

For the REW(Lya + N v)-log/? correlation, it appears that 
there are a few outliers with exceptionally weak Lya + N V 
lines (see the upper panel of Fig. [6j; they were identi- 
fied as weak-line quasars (WLQs) by Diamond-Stanic et al. 
(2009), including one of our Chandra Cycle 12 targets 
(J1412+0624). Another z > 4 HRLQ (J1659+2101) in our 
sample has REW(Lya + N v) values (17.9 A; see Table [HJl 
close to the WLQ criterion of REW(Lya + N v)< 15.4 A. 
These two objects also have weak C IV emission lines. Their 
REW(C IV) values are below the 3(7 negative deviation from 
the mean REW(C IV) of SDSS quasars (Wu et al. 2012). 
The weak emission lines of these objects may suggest they 
have a relativistic ally boosted optical/UV continuum which 
dilutes the emissio n line s. However, the broad-band SEDs of 
these objects (see £|4. 4i do not show the signature of a rel- 
ativistically boosted continuum (e.g., a parabolic profile of 
the SED; Nieppola et al. 2006). It is also possible that the 
weakness of their emission lines is caused by extreme quasar 
disk- wind properties as for the radio-quiet weak-line quasars 
(e.g., Richards et al. 2011; Wu et al. 2011, 2012). One of our 
moderately radio-loud quasars (J 1639+4340) also has a weak 
C IV line (see Table HQ. This object is more likely to fit in the 
latter scenario; and thus the apparent high NTD value of this 
object does not indicate that it has substantial jet-contributed 
optical/UV continuum. 

The jet-linked X-ray emission of RLQs is also likely 
beamed, although probably to a lesser extent compared to jet 
radio emission (e.g., see §6 of Mil). Therefore, the relation 
between emission-line REWs and relative X-ray brightness 
is also worth investigating. A correlation analysis was per- 
formed on a sample combining our RLQs and the objects in 
Mil that have available Lya + N V and C IV REW measure- 
ments (see Fig.|7]for the relation between line REWs and rel- 
ative X-ray brightness of this sample) using Spearman's rank- 
order analysis in the ASURV package; Spearman's rank-order 
analysis is usually preferred over Kendall's T test for samples 
where the number of objects is N ^ 300 The results are 
shown in Table [TOl Only a marginal correlation (93.3%) was 
found between REW(C IV) and AOox.rqq, which is consis- 
tent with the result in Richards et al. (201 1) that RLQs with 
weaker C IV lines also tend to be weaker in X-rays. However, 
the bottom-left panel (Aa ox ,RQQ vs. C iv REW) of Fig. [7] 
shows a large scatter in this correlation. No other significant 
correlations were found between relative X-ray brightness and 
emission-line strength for the RLQs (see TablefToland Fig. LZJ. 

4.3. X-ray Spectral Properties 

29 See Appendix A3. 3 of the ASURV manual at 
http://astrostatistics.psu.edu/statcodes/asurv. In our work, we choose 
Spearman's rank-order analysis when sample size N ^ 30 and Kendall's T 
test when N < 30. 
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Our Chandra Cycle 12 targets do not have sufficient counts 
for individual X-ray spectral analysis, and thus we obtained 
their effective power-law photon indices based on their band 
ratios (see Sj3]). The X-ray power-law photon index for RLQs 
is known to be anti-correlated with log/? (e.g., Wilkes & Elvis 
1987; Reeves & Turner 2000; Lopez et al. 2006; Saez et al. 
201 1). Our sample of HRLQs at z > 4 is generally consistent 
with this anti-correlation (see Fig. [HJ. Kendall's T test on a 
sample combining our HRLQs and the moderately radio-loud 
quasars at z > 4 in Saez et al. (201 1) shows a 2a-level anti- 
correlation between Tx and log/? (95.6% correlation proba- 
bility). The best-fit correlation for this sample, obtained via 
the EM linear-regression algorithm in the ASURV package.F1 
is 

r x = (1.95 + 0.28) + (-0.11 ±0.10)log/?. (8) 

This Tx-log/? relation (the dashed line in Fig. [8} is consistent 
with that for the z > 2 RLQs in Saez et al. (2011; see their 
Equation 1) and with that for the z < 2 RLQs in Reeves & 
Turner (2000; see Equation 2 of Saez et al. 2011), suggest- 
ing that this relation does not evolve with redshift. The two 
moderately radio-loud quasars observed in Chandra Cycle 12 
(J074 1+2520 and J 1639+4340), which do not satisfy our cri- 
terion for HRLQs, are outliers from this correlation. Instead 
of having the softer X-ray spectra generally associated with 
lower radio loudness, they have the hardest spectra among our 
objects. 

HRLQs often appear to have intrinsic X-ray absorption with 
Nh > 10 22 cm" 2 (e.g., Yuan et al. 2006). Bassett et al. (2004) 
found evidence of a similar level of intrinsic X-ray absorp- 
tion for z > 4 moderately radio-loud quasars. To assess pos- 
sible indications of intrinsic absorption for our Chandra Cy- 
cle 12 targets, we performed joint spectral analyses. We di- 
vided our six targets into two groups with log/? > 2.5 (four 
objects) and log/? < 2.5 (2 objects), respectively. For each ob- 
ject, the X-ray spectrum in observed-frame 0.5-8.0 keV (cor- 
responding to rest-frame w 2-40 keV) was extracted from a 
3"-radius circular region centered on the X-ray position using 
the standard CIAO routine PSEXTRACT. Background spec- 
tra were extracted from annular regions with inner and outer 
radii of 6" and 9", respectively. All background regions are 
free of X-ray sources. Joint X-ray spectral fitting was per- 
formed with XSPEC V12.6.0 (Arnaud 1996). We applied the 
C-statistic (Cash 1979) in the spectral fitting instead of the 
standard % 2 statistic because the C-statistic is well suited to 
limited X-ray count scenarios (e.g., Nousek & Shue 1989). 
We fit the spectra jointly using two models: (1) a power-law 
model with a Galactic absorption component represented by 
the wabs model (Morrison & McCammon 1983); (2) another 
model similar to the first, but adding an intrinsic (redshifted 
to the source rest frame) neutral absorption component, rep- 
resented by the zwabs model. Each source was assigned its 
own values of redshift and Galactic neutral hydrogen column 
density (see Column 4 of Table HJ. Table [Til shows the best- 
fit spectral parameters along with the errors or upper lim- 
its at the 90% confidence level for one parameter of interest 
(AC = 2.71; Avni 1976; Cash 1979). The best-fit power-law 
photon index values are consistent with those from band-ratio 
analysis (see Table |3). Our joint-fitting procedure is not able 
to provide tight constraints on the intrinsic X-ray absorption 
for either highly or moderately radio-loud quasars (see the up- 

30 This linear-regression algorithm is applicable to censored data. Our 
sample has two objects with upper limits upon Fx- 



per limits on intrinsic Nh in Table [TTJ mainly due to the lim- 
ited number of X-ray counts and their high redshifts. Adding 
an intrinsic absorption component does not improve the qual- 
ity of the joint fits. 

As discussed above, our two moderately radio-loud quasars 
appear to have a harder average X-ray spectrum than that 
of our HRLQs, which appears inconsistent with the known 
anti-correlation between X-ray power-law photon index and 
radio loudness. Possible heavy X-ray absorption may ex- 
ist in these two objects. We estimate the potential intrinsic 
X-ray absorption for these two moderately radio-loud quasars 
by assuming their unabsorbed power-law photon indices fol- 
low the log/?— Fx relation (see Equation 8; Fx, U nabs 

= 1.79 

for J0741+2520, r X;U nabs = 1-74 for J1639+4340), and then 
adding an intrinsic X-ray absorption component to make the 
apparent Tx agree with the values obtained from the band- 
ratio analysis (r x = 1.21 for J0741+2520, T x = 1.19 for 
J1639+4340; see Tabled. The required intrinsic hydrogen 
column densities are N H = 3.3 x 10 23 cm" 2 for J074 1+2520 
and N H = 2.0 x 10 23 cm" 2 for J1639+4340, which are about 
an order of magnitude higher than the values found for high- 
redshift RLQs in previous studies. The intrinsic X-ray ab- 
sorption in RLQs is often associated with absorption features 
in their optical/UV spectra (e.g., Elvis et al. 1998). The op- 
tical/UV spectrum of J1639+4340 (see the top-right panel of 
Fig. O shows several C IV and Si IV absorption features, al- 
though none of them has sufficient width to be identified as 
a broad absorption line (BAL). The optical/UV spectrum of 
J0741+2520 (see Fig. 3 of McGreer et al. 2009) does not 
show Si IV absorption features; it does not have coverage of 
the C IV region. 

4.4. The Radio -to -J-r ay Spectral Energy Distributions 

To investigate the broad-band SEDs of our RLQs, we gath- 
ered photometric data from the following sources: 

1. Radio: The 1.4 GHz flux densities are from the FIRST 
or NVSS surveys; the 5 GHz values are from the GB6 
or PMN catalogs or individual observations (see Sj3]); 
the flux densities at other frequencies were retrieved 
from the NED. 

2. Mid-infrared: From the all-sky catalog of the Wide-field 
Infrared Survey Explorer (WISE; Wright et al. 2010)0 
All of our objects have WISE detections in at least two 
bands (Wl at 3.4 /xm and W2 at 4.6 /im) except for 
GB 1428+4217 and PMN J 195 1+0 134. We examined 
the WISE images of these two objects. There appears to 
be a source at the optical location of GB 1428+4217, 
but it is strongly blended with a nearby brighter source 
9" away. There is no source at the optical location of 
PMN J1951+0134 in the relevant WISE images. All 
other objects are free of issues such as source blending, 
confusion, or mismatching. The high detection fraction 
(> 17/ 19) in the mid-infrared band is remarkable given 
the high redshifts of our objects. 

3. Near-infrared: From the Two Micron All Sky Survey 
(2MASS; Skrutskie et al. 2006). For the objects hav- 
ing SDSS spectroscopy, we also checked the SDSS 
DR7 quasar catalog which provides additional deeper 

31 See the WISE all-sky data release at 

http://wise2.ipac.caltech.edu/docs/release/allsky/ 
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2MASS photometry. Only our two moderately radio- 
loud quasars (J074 1+2520 and J 1639+4340) have 
2MASS detections, likely owing to their exceptional 
overall luminosities. GB 1428+4217 has additional 
near-infrared observations with the United Kingdom 
Infra-Red Telescope (UKIRT) reported by Fabian et al. 
(1999). 

4. Optical: From the SDSS database and/or the NED. The 
photometry for the bands that are seriously affected by 
the Lya forest is omitted. We do not include any UV 
photometric information for the same reason. 

5. X-ray: From this work. 

6. y-ray: The upper limits upon y-ray luminosities are 
derived from the two-year survey data of the Large 
Area Telescope (LAT; Atwood et al. 2009) onboard the 
Fermi Gamma-Ray Space Telescope. See below for de- 
tails. 

The radio-to- y-ray SEDs of our objects are plotted in Fig. [9] 
ordered by log/? in descending order. We compare the SEDs 
of our high-redshift HRLQs to those of similar RLQs at lower 
redshifts. There are two commonly used composite SEDs of 
RLQs in the literature: Elvis et al. (1994; E94 hereafter) and 
Shang et al. (2011; Sll hereafter). The E94 SED is biased 
toward X-ray bright quasars as they required their objects to 
have X-ray observations with good signal-to-noise ratio (see 
discussion in §3 of E94 and §5.1 of Richards et al. 2006a). 
For the Sll SED, although their sample is not X-ray selected, 
the majority of their RLQs have lower optical/UV luminosi- 
ties than those of our objects. The optical/UV luminosity can 
affect the shape of the broad-band SED (see Eqn. 5 for the de- 
pendence of otox on optical/UV luminosity). The a ox param- 
eter of RLQs also depends on radio luminosity (see Eqn. 5). 
Therefore we selected 10 RLQs in the Sll sample with com- 
parable optical/UV luminosity (logAZ^ [3000 A]> 45.9) and 
radio-loudness values (2.9 < log/? < 3.7) to those of ourz > 4 
HRLQs. We construct a comparison composite SED using 
these 10 objects (see the solid lines in Pig - [9j . We normalized 
the comparison SED to the observed data at rest-frame 2500 A 
for the following reasons: (1) rest-frame 2500 A is within 
a region without strong emission lines; (2) the optical/UV 
continuum generally does not have a substantial contribution 
from the jet (except per haps for a few objects with weak 
emission lines; see ^4.21 >; and (3) when SEDs are normal- 
ized at rest-frame 2500 A, the relative positions of X-ray data 
points directly reflect the a ox values of the relevant quasars. 
One caveat about the SEDs is that our multi-band photome- 
try is non-simultaneous, and some HRLQs are known to have 
strong variability. 

The 10 RLQs from Sll used in constructing our com- 
parison SED satisfy our log/? criterion for HRLQs; they 
are HRLQs at low redshifts (z < 1.4). The majority of 
our HRLQs (11/17) have higher X-ray luminosities than 
those of the comparison SED. Meanwhile, they do not show 
weaker optical/UV emission relative to their infrared emis- 
sion. Therefore, t he ge nerally higher Aa ox values of our 
z > 4 HRLQs (see §3jJ do reflect stronger X-ray emission 
instead of fainter optical/UV emission. This is an addi- 
tional illustration of the X-ray emission enhancement of z > 4 
HRLQs compared to HRLQs at lower redshift. The SEDs 
of RX J1028-0844, GB 1428+4217, and GB 1508+5714 



all reach maximal vLv values (over the 10 9 — 10 19 Hz range) 
in the X-ray band. This type of SED has been modeled as 
the emission generated by Compton scattering of synchrotron 
emission from the jet itself and/or of a powerful external pho- 
ton field such as the thermal emission of the accretion pro- 
cess (SSC/EC; see Fabian et al. 1998, 1999; Ghisellini et al. 
1998). The SED of Q 0906+6930 has also been described 
with a similar model (Romani 2006). Q 0906+6930, along 
with PMN J0525-3343 and GB 1713+2148, have notably 
stronger mid-infrared emission than that of the comparison 
SED, especially in the longer wavelength bands (W3 at 12 fim 
and W4 at 22 fim). This mid-infrared excess may have a con- 
tribution from the jet synchrotron radiation (e.g., Ghisellini 
et al. 1998). An alternative origin for the mid-infrared excess 
is the thermal emission of circumnuclear dust predicted by the 
clumpy torus model (e.g., Nenkova et al. 2008). The SEDs 
of our HRLQs are significantly different from those of typi- 
cal BL Lac objects which generally have a parabolic profile 
peaking in the near-infrared band (e.g., Nieppola et al. 2006). 

Many blazars are high-energy y-ray emitters above 
100 MeV (e.g., Hartman et al. 1999). Q 0906+6930 has 
a tentative claimed detection (~ 1.5(7) in its y-ray obser- 
vation by the Energetic Gamma-Ray Experiment Telescope 
(EGRET) onboard the Compton Gamma-Ray Observatory 
(Romani et al. 2004). The Fermi LAT has an order of mag- 
nitude better sensitivity and positional accuracy than EGRET. 
Romani (2006) suggested Q 0906+6930 should be detected 
by the Fermi LAT in its first year of operation. However, al- 
though the first two-year survey of the Fermi LAT has de- 
tected > 800 blazars at z < 3 (e.g., Ackermann et al. 2011), 
none of our HRLQs, including Q 0906+6930, is detected at 
present. We searched the second source catalog of the Fermi 
LAT (2FGL; Nolan et al. 2012) and found no counterparts for 
our HRLQs within a 10' matching radius To obtain upper 
limits upon y-ray luminosity, we first estimated the photon 
flux limits (in units of photons cm~ 2 s _1 ) for all objects based 
on the point-source flux limit map of the Fermi LAT two-year 
survey (see Fig. 1 of Ackermann et al. 2011). We converted 
the photon flux limits to energy flux limits between observed- 
frame 100 MeV-100 GeV and then flux-density limits at rest- 
frame 1 GeV, assuming a typical power-law photon index 
of FermZ-detected FSRQs (r^ = 2.4; see Fig. 18 of Acker- 
mann et al. 2011). The upper limits upon the monochro- 
matic y-ray luminosity at rest-frame 1 GeV, calculated from 
the flux-density limits at the same frequency, are in the range 
of(1.2-2.6)xl0 47 ergs- 1 Hz" 1 (see Tableland Fig.©. For 
some of our z > 4 HRLQs with strong X-ray enhancements 
(e.g., Q 0906+6930, RX J1028-0844, GB 1428+4217, and 
GB 1508+5714), the y-ray luminosity limits appear inconsis- 
tent with the prediction of the synchrotron-Compton blazar 
emission model in Romani et al. (2006; see their Fig. 3). Fur- 
ther Fermi LAT survey observations will provide tighter y-ray 
luminosity upper limits (by a factor of ~ 2) or perhaps detec- 
tions of z > 4 blazars, which will help constrain the emission 
models for blazars at high redshift. Sbarrato et al. (2012a) 
found a positive correlation between the luminosity of op- 
tical/UV broad emission lines and the y-ray luminosity for 
ferazi-detected blazars. This correlation may suggest that the 
z > 4 blazars with higher optical/UV emission-line luminos- 
ity should be easier to detect with Fermi. 

32 We also searched for the z > 4 HRLQs without sensitive X-ray coverage 
listed in Table [2] and found none of them having counterparts in the 2FGL 
catalog. 
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Six objects with SDSS spectroscopy have available black- 
hole mass (Mbh) and Eddington-ratio (L/Ledd) estimates (see 
Table [H in the catalog of Shen et al. (2011)0 The bolo- 
metric luminosity in Shen et al. (2011) was calculated from 
the optical/UV luminosity and the bolometric correction from 
the composite SED of Richards et al. (2006a). Therefore, 
the bolometric luminosity and the Eddington ratio only ac- 
count for the emission of the quasar from its accretion activ- 
ity (i.e., excluding beamed jet emission). All of the six ob- 
jects have M BH > 1O 9 M except for GB 1508+5714. These 
five objects have estimated Eddington ratios between 0.1 and 
1, which are typical for SDSS quasars at high redshift (e.g., 
Shen et al. 2008). GB 1508+5714 has a smaller M BH but 
an Eddington ratio L/Lgdd ~ 3. This object is similar to the 
Ferm/-detected FSRQs in Shaw et al. (2012) which have 
smaller Mbh and higher Eddington ratios than those of opti- 
cally selected quasars. The objects in Shaw et al. (2012) show 
significant non-thermal (synchrotron) emission in the optical 
band associated with the jet. However, GB 150 8+57 14 has a 
smaller NTD value (consistent with unity; see ^4. 2\ than the 
majority of the objects in Shaw et al. (2012; see their Fig. 3), 
which suggests that this object does not have substantial jet- 
linked optical/UV continuum emission. It is worth noting that 
the Mbh measurements for our objects are all based on the 
C IV emission line due to their high redshifts. The C IV esti- 
mator for Mbh is more affected by the quasar disk wind than 
the H/3 and Mg II estimators, and this may lead to errors for 
individual Mbh measurements (e.g., Shen et al. 2008). Thus, 
the apparent super-Eddington emission of GB 1508+5714 
should not be over-interpreted. 

5. SUMMARY AND FUTURE STUDIES 

We have compiled a sample of 17 highly radio-loud quasars 
(log/? > 2.5) at z > 4 with sensitive X-ray coverage by Chan- 
dra, XMM-Newton, or Swift, Four of them were targeted in 
new Chandra Cycle 12 observations, while the other 13 ob- 
jects have sensitive archival X-ray coverage. Our sample of 
HRLQs represents the top ~ 5% of the total RLQ popula- 
tion in terms of radio loudness. They are among the most- 
luminous quasars in both the radio and optical/UV bands. 
The X-ray and broad-band properties of our HRLQs are pre- 
sented and investigated. We have also reported new Chandra 
Cycle 12 observations of two moderately radio-loud quasars 
(1 < log/? < 2.5) at z > 4 which are among the most optically 
luminous RLQs. Our main results are summarized as follows. 

1. All of our Chandra Cycle 12 targets are detected 
in X-rays. None of them shows detectable ex- 
tended X-ray emission. All of the archival ob- 
jects in our sample of z > 4 HRLQs were also de- 
tected by Chandra, XMM-Newton, or Swift, includ- 
ing four objects for which their archival Chandra 
(GB 1713+2148),XMM-Afewfon(SDSS J1235-0003), 
or Swift (PMN Jl 155-3107 and PMN J1951+0134) 
observations are first reported in this work. See $3] 

2. Our HRLQs at z > 4 show stronger X-ray emission (by 
a typical factor of w 3) than HRLQs at lower redshift 

33 Another HRLQ in our sample, J 1026+2542, has M B h estimated by 
Sbarrato et al. (2012b) based on the peak frequency of disk emission and 
the total disk luminosity. The estimated range is M B H=(1.8^t.5)xl0 9 M 1: ;. 
Although this object has SDSS spectroscopy, its Mbh cannot be measured 
based on available emission-line widths because none of the standard lines 
(H/3, Mg II, and C IV) is covered due to its high redshift (z = 5.304). 



with similar optical/UV and radio luminosities. This 
contrasts with the behavior of moderately radio-loud 
quasars at z > 4, implying that the high-redshift X-ray 
enhancement occurs within the jet-linked component 
that becomes increasingly prominent for HRLQs. We 
examined possible biases in our analysis and found that 
our result remains robust after controlling for a variety 
of potential issues. A similar X-ray emission enhance- 
ment is also found for z = 3—4 HRLQs, which provides 
statistically independent evidence for our findings. The 
observed X-ray enhancement is not likely to be caused 
by the X-ray spectral cur vature found in some high- 
redshift HRLQs. See gjj 

3. Our z > 4 HRLQs are generally consistent with the 
positive correlations between optical/UV emission-line 
REWs (Lya + N V and C IV) and radio loudness found 
for the typical RLQ population. These positive correla- 
tions suggest that the optical/UV continua of HRLQs 
usually have a negligible contribution from the rela- 
tivistic jets (except for perhaps a few outliers with weak 
emission lines). Our z > 4 HRLQs do not show strong 
correlations between the relative X-ra y bri ghtness and 
optical/UV emission-line REWs. See S 14.2I 

4. Our sample of z > 4 HRLQs generally follows the 
known anti-correlation between X-ray power-law pho- 
ton index and radio loudness. However, the two tar- 
geted moderately radio-loud quasars, having the hard- 
est X-ray spectra among our objects, appear to be out- 
liers from this correlation, suggesting the possible pres- 
ence of intrinsic X-ray absorption (Nn ~ 10 23 cm -2 ). 
See g3] 

5. The majority of z > 4 HRLQs have higher X-ray lumi- 
nosities than those of the matched low-redshift HRLQs 
used to construct a comparison SED, which further il- 
lustrates and supports the X-ray emission enhancement 
of z > 4 HRLQs over those at lower redshifts. The 
SEDs of our objects with the highest X-ray luminosities 
are consistent with the expectations from SSC/EC mod- 
els. Some of our HRLQs also show an excess of mid- 
infrared emission which may be contributed by syn- 
chrotron radiation from the relativistic jets. None of 
our HRLQs has been dete cted by the Fermi LAT in its 
two-year survey. See ^4.41 

We have performed the first systematic X-ray survey of 
HRLQs at z > 4. Our primary new result is the measurement 
of an increase (by a factor of ss 3) in the relative X-ray emis- 
sion of z > 4 HRLQs compared to their low-redshift counter- 
parts. This work adds the relative X-ray brightness of HRLQs 
to the RLQ properties that evolve with redshift, such as the 
frequency of intrinsic X-ray absorption (see §1) and perhaps 
the RLQ fraction (e.g., Jiang et al. 2007; Singal et al. 2012). 
The number density of hard X-ray selected FSRQs has been 
found to evolve strongly with redshift, with a peak number 
density of Zpeak ~ 

3-4 (Ajello et al. 2009). Unfortunately, 
the sample utilized by Ajello et al. (2009) from the Swift 
Burst Alert Telescope (BAT; Barthelmy et al. 2005) does not 
have any FSRQ detections at z > 4, and thus there is consid- 
erable uncertainty in making quantitative comparisons at the 
high redshifts most applicable to our work. We do note, how- 
ever, that the SED evolution we find for high-redshift HRLQs 
(which are primarily FSRQs) will affect the interpretation of 
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number-density evolution results for X-ray selected samples. 
For example, it may help to explain the higher z pea k value 
found for FSRQs than for general quasars or X-ray selected 
AGNs. 

While the limited increase in X-ray luminosities confirms 
earlier findings that the small-scale X-ray jet-linked emission 
is not dominated by IC/CMB emission, our results are consis- 
tent with a fractional contribution from IC/CMB (~ 6% of the 
nuclear X-ray emission at z = 1.3) rising with redshift. This 
does not conflict with growing evidence that RLQ jets have 
mildly relativistic bulk velocities beyond several kpc, as in- 
ferred from radio core/jet prominences (e.g., yj e t ~ 1.2 — 1.5; 
Mullin & Hardcastle 2009), or that various predictions of 
IC/CMB models are not obviously met in X-ray observations 
of large-scale RLQ jets (e.g., Hardcastle 2006; Marshall et 
al. 2011; Massaro et al. 2011; note these studies are based 
almost exclusively on RLQs with z < 2.5). It may be that 
IC/CMB emission is relevant only over a limited range of spa- 
tial scales (e.g., w 1-5 kpc), with quasar-related photon fields 
dominating at smaller distances (e.g., Ghisellini & Tavecchio 
2009) and the jet decelerating at larger distances. For exam- 
ple, a drop in the bulk jet velocity from j8 = 0.995 to j3 = 0.75, 
corresponding to a decrease from y = 10 at < 5 kpc to y = 1.5 
at large scales, would act to diminish X-ray IC/CMB emis- 
sion by a factor of w 44. Jet deceleration appears to be re- 
quired in at least some RLQs even in the context of IC/CMB 
modeling of the large-scale X-ray jet emission, when coupled 
with VLBI constraints (Hogan et al. 2011). Such decelera- 
tion could be partially due to mass entrainment; simulations 
predict that even in powerful jets the bulk Lorentz factor de- 
creases by a factor of ~2 on kpc scales (Bowman et al. 1996). 
Although faster jets are more stable against Kelvin-Helmholtz 
instabilities (e.g., Perucho et al. 2005), the gradual accumula- 
tion of velocity-shear perturbations may similarly act to slow 
the jet. The Compton-rocket effect (e.g., Ghisellini & Tavec- 
chio 2010) provides another means of draining jets. Whatever 
the dominant mechanism, the postulated degree of deceler- 
ation could plausibly be accomplished without disrupting or 
decollimating the jet; for example, in their classic work Kaiser 
& Alexander (1997) estimate that a powerful jet with y e t = 2 
is persistently stable (although see also Wang et al. 2011b). 
This deceleration possibility can be tested with additional ob- 
servations of large-scale jets in high-redshift RLQs0 In a 
fractional IC/CMB scenario, the enhanced X-ray emission 
from z > 4 HRLQs is solely due to the increased energy den- 
sity of the CMB; the jet bulk velocities, the disk/corona ac- 
cretion structure, and the host-galaxy properties of individual 
RLQs are not required to evolve with redshift. 



Alternative interpretations of the observed X-ray emission 
enhancement of z > 4 HRLQs are also possible. For example, 
increased host-galaxy star formation could provide additional 
seed photons for external Compton upscattering in the jet (see 
§4.1.3). Also, a decrease in the jet bulk Lorentz factor at 
high redshift (a possibility suggested by Volonteri et al. 201 1) 
might increase the ratio of X-ray to radio emission, particu- 
larly in a multi-component jet. It does not appear likely that 
X-ray spectral curvature drives our results, as the observed- 
frame ultra-soft luminosities are also enhanced. After careful 
consideration of selection effects, there is no obvious indica- 
tion that our sample contains a disproportionate number of 
"extreme" HRLQs, with properties (low inclinations, strong 
variability) likely to correlate with enhanced X-ray emission. 

Future X-ray observations of additional z > 4 HRLQs, e.g., 
a Chandra snapshot survey of the 11 objects listed in Ta- 
ble 2, would further improve the sample statistics and bet- 
ter constrain the possible mechanisms for the X-ray enhance- 
ment. For example, the HRLQs at z = 3-4 apparently have a 
similar level of X-ray enhancement to that of z > 4 HRLQs, 
which would not agree with the expectations from the frac- 
tional IC/CMB model. However, more objects are needed 
to reduce the large uncertainties on the mean Aa ox values in 
current analyses. Identifying potential intrinsic X-ray absorp- 
tion in our objects, especially for the two moderately radio- 
loud quasars (J0741 +2520 and J 1639+4340), requires deeper 
X-ray spectroscopic observations. The ongoing Fermi LAT 
survey will provide tighter y-ray luminosity upper limits or 
perhaps y-ray detections of z > 4 HRLQs, which will be use- 
ful for modeling their SEDs and for constraining the relative 
strengths of the jet emissions from different mechanisms. 
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TABLE 1 
X-ray Observation Log 



Object Name 


R.A. (J2000) 
(deg) 


Dec (J2000) 
(deg) 




Detector 


Observation 
Date (UT) 


Observation 
ID 


Exp. Time b 
(ks) 


Reference 


Gxu.Yid.Yii Cycle 12 Objects 


















SDSS J074154.71 + 252029.6 


115.4780 


25.3416 


5.194 


ACIS-S 


2010 Oct 04 


12171 


4.0 


1,2 


SDSS J102623.61 +254259.5 


156.5984 


25.7165 


5.304 


ACIS-S 


2011 Mar 10 


12167 


5.0 


1 


SDSS J141209.96 + 062406.9 


213.0415 


6.4019 


4.467 


ACIS-S 


2011 Mar 16 


12169 


4.0 


1 


SDSS J142048.01 + 120545.9 


215.2000 


12.0962 


4.027 


ACIS-S 


2011 Mar 20 


12168 


4.0 


1 


SDSS J163950.52 + 434003.6 


249.9605 


43.6677 


3.980 


ACIS-S 


2010 Aug 08 


12170 


4.0 


1 


SDSS J165913.23 + 210115.8 


254.8051 


21.0211 


4.784 


ACIS-S 


2011 Jan 22 


12172 


6.5 


1 


Archival X-ray Data Objects 


















rSo Ulzl + U34 / 


in icon 


3.7851 


A 1 in 
4.1 JU 


ACIS-S 


zUUZ reb U/ 


T 1 CI 

JlDl 


5.7 


1 
J 


PMN TnnA 9Q1S 
rlvlrN JUjZ^t — zylo 


c. 1 1 








9nrn Hph i f\ 

ZUUj ucl io 




J.O 


-T.J 


PMN J0525 - 3343 e 


81.2758 


-33.7182 


4.401 


MOS f 


2001 Sep 15 


0050150301 


19.0 


4,6,7 


Q 0906 + 6930 


136.6281 


69.5085 


5.480 


ACIS-S 


2005 Jul 1 


5637 


30.0 


8.9 


SDSS J091316.55 + 591921.6 d,g 


138.3190 


59.3227 


5.122 


ACIS-S 


2002 Mar 7 


3034 


9.8 


10 


RXJ1028.6-0844 


157.1617 


-8.7441 


4.276 


MOS f 


2003 Jun 13 


0153290101 


21.3 


11,12,13 


PMN Jl 155 -3107 


178.7632 


-31.1330 


4.300 


XRT 


2008 Jan 07 


00036791002 


4.2 


1,14 


SDSS JI23503.03 - 000331.7''- 


188.7626 


-0.0588 


4.673 


MOS 


2010 Jul 01 


0651740301 


5.3 


1 


CLASS J1325 + 1123 d - E 


201.3021 


11.3916 


4.415 


ACIS-S 


2003 Mar 02 


3565 


4.7 


15 


GB 1428 + 4217 dc 


217.5989 


42.0768 


4.715 


MOS f 


2003 Jan 17 


0111260701 


14.2 


16,17,18,19 


GB 1508 + 5714 d >6 


227.5122 


57.0454 


4.313 


ACIS-S 


2001 Jun 10 


2241 


90.1 


20,21,22,23 


GB 1713 + 2148 d 


258.8385 


21.7588 


4.011 


ACIS-S 


2004 Jun 07 


4815 


9.5 


1,3,16 


PMN J1951 + 0134 


297.9001 


1.5785 


4.114 


XRT 


2007 Mar 15 


00036263002 


7.8 


1,14 


a Redshift for each source. The redshift values for objects in the SDSS DR7 quasar catalog are from Hewett & Wild (2010). 

The exposure times are corrected for detector dead time. 
c References: (1) This work; (2) McGreer et al. (2009); (3) Vignali et al. (2003a); (4) Hook et al. (2002); (5) Lopez et al. (2006); (6) Fabian et al. (2001); (7) Worsley et al. (2004a); (8) 
Romani et al. (2004); (9) Romani (2006); (10) Vignali et al. (2003b); (1 1) Zickgraf et al. (1997); (12) Grupe et al. (2004); (13) Yuan et al. (2005); (14) Healey et al. (2008); (15) Bassett 
et al. (2004); (16) Hook & McMahon (1998); (17) Fabian et al. (1998); (18) Fabian et al. (1999); (19) Worsley et al. (2004b); (20) Siemiginowska et al. (2003); (21) Yuan et al. (2003); (22) 
Cheung (2004); (23) Yuan et al. (2006). 
d These archival sources are covered by the SDSS footprint. 

e These sources were observed in multiple XMM-Newton or Swift observations. The parameters listed in this table are for the observation with the longest MOS or XRT exposure time. MOS 



exposure times have been corrected for the background flares. 

^ These objects were observed by both the MOS and pn detectors. We list the MOS exposure times here. 
s These archival objects have SDSS spectroscopy; they appear in the SDSS DR7 quasar catalog. 



TABLE 2 

HRLQS AT z > 4 WITHOUT SENSITIVE X-RAY COVERAGE 



Object Name 


R.A. (J2000) 


Dec (J2000) 




m, b 


Mi 


/l.4 GHz C 


logR 




(deg) 


(deg) 








(mJy) 




SDSS J003126.79+ 150739.5 


7.8616 


15.1276 


4.297 


20.0 


-27.2 


41.0 


2.824 


SDSS J030437.21 +004653.5 


46.1550 


0.7815 


4.260 


20.2 


-27.1 


20.0 


3.259 


SDSS J081333.32 + 350810.8 


123.3888 


35.1363 


4.945 


19.2 


-28.2 


20.0 


2.663 


NVSS J10501 1-044254 


162.5480 


-4.7150 


4.270 


23.4 


-23.8 


6.0 


3.338 


NVSS Jl 12310 -215405 


170.7921 


-21.9016 


4.110 


22.9 


-24.1 


48.0 


4.041 


SDSS J123142.17 + 381658.9 


187.9257 


38.2830 


4.128 


20.1 


-26.9 


25.0 


2.777 


SDSS J123726.26 + 65 1724.4 


189.3594 


65.2901 


4.300 


20.5 


-26.6 


24.0 


2.785 


SDSS J124230.58 + 542257.3 


190.6274 


54.3826 


4.766 


19.7 


-27.7 


19.0 


2.678 


PMN J2134- 0419 


323.5501 


-4.3194 


4.350 


20.1 


-27.1 


333.0 


3.714 


SDSS J222032.50 + 002537.5 


335.1354 


0.4271 


4.218 


19.9 


-27.1 


61.0 


3.531 


PMN J2314 + 0201 


348.7030 


2.0309 


4.110 


20.2 


-26.9 


117.0 


3.308 



il The redshift values for SDSS objects, all of which appear in the SDSS DR7 quasar catalog, are from Hewett 
& Wild (2010). 

b The apparent i'-band magnitude. For objects not in the SDSS footprint, the (-band magnitude was estimated 
from their fi-band magnitude by assuming a power law with a spectral index of a v = —0.5. 
c Radio flux density at observed-frame 1.4 GHz from FIRST or NVSS. 
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TABLE 3 

X-ray Counts for Objects without Previously Published X-ray Photometry 





Full Band 


Soft Band 


Hard Band 


Band 




Object Name 


(0.5-8.0 keV) a 


(0.5-2.0 keV) a 


(2.0-8.0 keV) a 


Ratio b 





Chandra Cycle 12 Objects 



SDSS J074154.71- 


-252029.6 


29.8; 


-6. 
-5. 


SDSS J102623.61 - 


-254259.5 


59.3; 


-8. 
-7. 


SDSS J141209.96 - 


- 062406.9 


29.9; 


6. 
-5. 


SDSS J142048.01 - 


- 120545.9 


48.9; 


1-8. 
-7. 


SDSS J163950.52- 


-434003.6 


42.9; 


(-7. 
6. 


SDSS J165913.23- 


-210115.8 


43.6; 


r7. 

6. 



17.83! 11.1'H 1.2118$ 

42.6+™ 14.3+« 0.34+0;" 1.73+°J? 

m.9tl A 3 9.9i« 0.52+°|7 1.33+°« 

33 9+ 6 ' 13 1+ 4 - 7 39+° 16 161+° " 

-"■ v -5.8 iJa -3.6 U - Jy -0.12 1 - D1 -0.32 

25.8+j! 15.7+5+ 061+ o.25 U9 +o.34 

30.8+|;f 11.1+11 °- 36 -o:l2 lJ5 -of 6 



Archival X-ray Data Objects 



PMNJ1 155 -3107 12.1+lSj 10.1+1] < 8.3 < 0.82 > 1.33 

SDSS J123503.03- 000331.7 25.9+g ! < 18.3 < 15.9 

GB 1713 + 2148 92.5+^j 7 63.2+?;° 25.9+f g 0.41+°;^ 1.63±g;g 

PMNJ1951+0134 45.2+™ 25.4+j : J 20.7+|f 0.81+° : ^ 1.54±jj'^ 



a Errors on the X-ray counts were calculated using Poisson statistics corresponding to the 1 o~ significance level according 
to Tables 1 and 2 of Gehrels (1986). 

b The band ratio is defined here as the number of hard-band counts divided by the number of soft-band counts. The 
errors on the band ratio correspond to the lo~ significance level and were calculated using equation (1.31) in §1.7.3 of 
Lyons (1991). The band ratios for all of the Chandra objects observed in the same cycle can be directly compared with 
one another. 

c The effective power-law photon indices were calculated using the PIMMS tool (version 3.9&). The effects of the 
quantum-efficiency decay over time at low energies of the ACIS detector were corrected for Chandra observed objects. 
The Chandra ACIS Cycle 5 response was used for GB 1713+2148, while the Chandra Cycle 12 response was used for 
the other objects. For Jl 155— 3107 and J1235— 0003 which do not have photon index estimation, we adopt a typical 
RLQ photon index value Y = 1.6 in the following analyses. 



TABLE 4 

X-RAY, OPTICAL/UV, AND RADIO PROPERTIES 



Object Name 
(1) 



(2) 



(3) 



Na 
(4) 



Count 
Rate c 
(5) 



Fx" 
(6) 



fl kcV C 
(V) 



(2 - 10 keV) f 
(8) 



(9) 



■^2500 A 



log L uv 
(2500 A)' 



log L r 
(5 GHz) k 



logfi 1 



AOOX.RQQ (0") m Aa ox , RLQ n 



"(10) 


(ID 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


1.95 


31.97 




33.16 


1.059 


-1.71 


0.07 (0.5) 


-0.02 


0.62 


31.49 


-0.38 


35.15 


3.542 


-1.31 


0.39 (2.7) 


0.07 


0.65 


31.40 




34.36 


2.698 


-1.51 


0.18 (1.2) 




0.56 


31.26 


-0.36 


34.55 


3.052 


-1.34 


0.33 (2.3) 




4.29 


32.14 


-0.32 


33.93 


1.621 


-1.80 


-0.01 (0.1) 


o.|j 


0.60 


31.40 




34.25 


2.557 


-1.39 


0.30 (2.1) 





Chandra Cycle 12 Objects 
SDSS J074154.71 + 252029.6 
SDSS J102623.61 + 254259.5 
SDSS J141209.96 + 062406.9 
SDSS J142048.01 + 120545.9 
SDSS J163950.52 + 434003.6 
SDSS J165913.23 + 210115.8 

Archival X-ray Data Objects 



18.54 


-29.04 


4 


26 


4.44 


20.03 


-27.50 


1 


80 


8.55 


19.44 


-27.74 


2 


1 1 


4.73 


19.80 


-27.18 


1 


72 


8.47 


17.78 


-29.13 


1 


36 


6.45 


20.26 


-27.17 


5 


47 


4.77 



PSS 0121+0347 




18 


57 


-28 


44 


3.19 


PMNJ0324-2918 




18 


65 


-28 


61 


1.19 


PMN J0525 - 3343 




18 


63 


-28 


52 


2.19 


Q 0906 + 6930 




19 


85 


-27 


76 


3.64 


SDSS J091316.55 + 


591921.6 


20 


39 


-27 


03 


3.85 


RX J1028.6-0844 




19 


14 


-27 


95 


4.60 


PMN Jl 155 -3107 




19 


28 


-27 


90 


6.04 


SDSS J123503.03 - 


000331.7 


20 


10 


-27 


20 


1.90 


CLASS J1325+ 1123 


19 


18 


-28 


01 


1.93 


GB 1428 + 4217 




19 


10 


-28 


18 


1.39 


GB 1508 + 5714 




19 


92 


-27 


16 


1.46 


GB 1713 + 2148 




21 


42 


-25 


53 


5.05 


PMN J1951 +0134 




19 


69 


-27 


40 


16.42 



-1.04 
+ 1.52 
30 
.36 
-1.08 
-.72 
-1.45 
.54 
26 
-.02 
-0.85 



10.30l| Io 
13.70+?$ 



10.72 
0.51 



-0.60 
+0.34 
-0.29 



5.04 



+ 1.26 
1.05 



+0.70 



42.06+ 069 
6 63+ ' 94 

aoJ -0.83 



1.89 
3.51 
1.93 
3.45 
2.57 
2.15 



6.80 
24.03 

7.73 
17.84 

8.67 
14.12 



45.20 
45.76 
45.15 
45.44 
45.13 
45.34 



, -,,+0.44 
1 - zl -0.37 
1 73+°- 32 
, ,1+0.43 
1 - 3J -0.38 
1 61 + 035 



1 iq+0.34 
, 7^+0.40 



3.19 


26.78 


45.64 


2 10+ 040 
z - lu -0.30 


1.77 


31.78 


-0.33 


34.47 


2.573 


-1.47 


0.28 (1.9) 


5.83 


39.70 


45.89 


LMJ -0.40 


1.73 


31.84 


+0.30 


34.92 


2.953 


-1.40 


0.35 (2.4) 


27.70 


159.40 


46.46 


1 67+ 002 

1 - D '-0.02 


1.73 


31.81 


+0.06 


34.84 


2.903 


-1.17 


0.58 (4.0) 


4.16 


24.78 


45.80 


1 6+° 1 
iD -0.1 


0.61 


31.51 


+0.17 


34.64 


3.008 


-1.31 


0.40 (2.7) 


0.16 


0.95 


44.34 


>0.71 


0.36 


31.23 


-0.67 


34.07 


2.717 


-1.76 


-0.09 (0.6) 


38.14 


179.40 


46.49 


1 40+° 03 
1 - w -0.03 


1.26 


31.65 


-0.30 


35.03 


3.329 


-1.09 


0.63 (4.3) 


5.63 


29.79 


45.34 


> 1.33 


1.01 


31.56 


+0.53 


34.37 


2.727 


-1.36 


0.36 (2.4) 


2.15 


10.57 


44.96 




0.16 


30.81 




33.87 


3.051 


-1.22 


0.39 (1.9) 


1.62 


10.71 


45.29 


, o n +0.50 


1.03 


31.59 


-0.09 


34.43 


2.722 


-1.53 


0.19 (1.3) 


70.07 


447.30 


46.95 


, 71+0.03 
1 ''-'-0.03 


1.15 


31.68 


+0.37 


34.87 


3.065 


-0.93 


0.80 (5.5) 


15.08 


76.68 


46.13 


1 55+ 005 

1 ' JJ -0.05 


0.24 


30.94 


+0.13 


34.93 


3.871 


-0.96 


0.67 (3.4) 


2.36 


12.40 


45.29 


1 63+ 024 


0.13 


30.61 


-0.30 


35.11 


4.501 


-1.16 


0.42 (2.1) 


8.78 


42.79 


45.51 


1 54+ 033 

1 ' J -0.30 


0.69 


31.36 


+0.24 


34.48 


3.036 


-1.23 


0.45 (3.1) 



0.033 
0.0© 
0.3th 

-0.5J 
-o.te 

0.3> 
0.2^ 



The apparent (-band magnitude. 

The absolute /-band magnitude, corrected for Galaelie extinction. 



The count rale of the ChaiiJnt-ob served sources in the observed- frame soft X-ray band (0.5 — 2.0 keV), 

d ^_ , , * — o c , m -14 „ rn[ „ 



The Galactic absorption-coiTected observed-frame flux between 0.5 — 2.0 keV in i 

6 The flux density at rest-frame 2 keV, in units of 10~^~ ergs cm - ' s — * Hz - * . 
f 

The logarithm ol the X-ray luminosity in the rest-lrame 2-10 keV band, eoneeted for Galai 
^ The X-ray power-law photon index (also see Tabled. 

^ The flux density at rest-frame 2500 A in units of 10~ 27 ergs em -2 s -1 Hz -1 . 

The logarithm of the monochromatic UV luminosity at rest-frame 2500 A. 
^ The radio spectral index a between observed -frame 1 .4—5 GHz, defined as /y « v a . 



n units oriO - 
-2 c -l 



The logarithm of monochromatic radio luminosity at rest-frame 5 GHz. 



'I lie logarithm ol" radii i loudness: see ; : 3 fur deli nil ion. 



Aa ox,RQQ 
i - 0.003-0.07. 
' Ao bx,RLQ : thedifft 



the difference between the measured Cf () x ihe expected U x for RQQs with 

Iween the measured Oj\ and the expected Ctox for RLQs with similar UV and rad: 



UV luminosity, defined by the ctox — ^-igrjQ ^ relation in equation (3) of Just ct ;U. (2007). The statistical significance of this difference, i 
isities, defined by the L 2 keV - ^2500 A ~~ L $ GHz relation in Table 7 of Miller et al. (2011) for their full sample. 



s measured in units of the RMS Oqx defined in Table 5 of Sleffen el al. (2006). The error range of Au x 
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TABLE 5 
Results of Peto-Prentice Tests 



Parameter 


iNO. 01 


No. of Comparison Objects 


Mean 


Mean 


oiausuc 


Null-hypothesis 




Our HRLQs 


(No. of X-ray Limits) 


(Our HRLQs) 


(Comparison Objects) 




Probability 


ACKox.RQQ 


12 


283(14) 


0.435 ±0.054 


0.251 ±0.008 


3.844 


1.21 x 10" 4 


ACKox.RLQ 


12 


283(14) 


0.1 67 ±0.049 


0.005 ±0.007 


3.814 


1.37 x ur 4 


Aoiox.RQQ (Worrall et al.'s Objects) 


8" 


81(7) 


0.475 ±0.065 


0.302±0.018 


2.682 


7.32 x 10~ 3 


Aa x,RLQ (Worrall et al.'s Objects) 


8 


81(7) 


0.204 ±0.059 


0.014 ±0.016 


3.181 


1.47 x ur 3 


Aa ox ,RQQ (FIRST-selected) 


8 


70(0) 


0.475 ±0.065 


0.241 ±0.011 


3.892 


9.94 x 10 -5 


Aa OXjRL0 (FIRST-selected) 


8 


70(0) 


0.205 ±0.059 


0.002 ±0.011 


3.814 


1.37 x 10~ 4 


AOox,RQQ (logi 5 GHz > 34) b 


12 


133(7) 


0.435 ±0.054 


0.294 ±0.012 


2.668 


7.63 x 10~ 3 


Aa ox ,RLO QogL 5 ghz > 34) 


12 


133(7) 


0.1 67 ±0.049 


-0.001 ±0.011 


3.628 


2.86 x Kr 4 


Aa x,RQQ ( lo g L 2500 A > 30 - 9)b 


12 


90(6) 


0.435 ±0.054 


0.305 ±0.016 


2.363 


1.81 x ur 2 


Aa ox ,RLQ ( lo g L 2500 A > 30 ' 9 ) 


12 


90(6) 


0.167 ±0.049 


0.018±0.013 


3.126 


1.77 x ur 3 


AOox,RQQ (a r > -0.5) 


U c 


90(2) 


0.458 ±0.054 


0.261 ±0.014 


3.782 


1.56 x 10~ 4 


Aa ox ,RLQ (Or > -0.5) 


11 


90(2) 


0.188 ±0.049 


0.029 ±0.013 


3.587 


3.35 x 10~ 4 



NOTE. — For the detailed definition of the Peto-Prentice test statistic see Feigelson & Nelson (1985). The null-hypothesis probability was calculated from the test statistic 
based on a Gaussian distribution, e.g., 1 — Pq — 1.21 x 10~ 4 , where Pq is the cumulative Gaussian probability at 3.844(7. 

a The number of radio-selected objects in our z > 4 HRLQ sample that satisfy m; < 20. One object, GB 1713+2148, has > 20, and thus is excluded from the two-sample tests, 
k The range of radio luminosity L5 ghz > 34 and optical/UV luminosity L s > 30.9 were determined based on our z > 4 HRLQs with m, < 20. Thus all 12 of our HRLQs 
surviving this m* cut are included in these tests. 

L All of our 12 z > 4 HRLQs with m; < 20 are FSRQs except for J1412+0624 which does not have available CK r measurement. 



TABLE 6 

Results of Peto-Prentice Tests for Additional Redshift Ranges 



Sample 


No. of HRLQs 


Mean(AOox.RQQ) 


Mean(Aa 0J[ .RLQ) 




(No. of X-ray Limits) 




3 sj z < 4 


9(1) 


0.468 ±0.056 


0.167 ±0.046 


Z < 3 


274(13) 


0.244 ±0.008 


0.000 ±0.007 


Test Statistic 




4.954 


5.428 


Null-hypothesis Probability 




7.27 x 10~ 7 


5.70 x ur 8 


2 < z < 3 


39(4) 


0.249 ±0.023 


-0.030±0.019 


z < 2 


235(9) 


0.243 ± 0.009 


0.005 ± 0.008 


Test Statistic 




0.176 


1.981 


Null-hypothesis Probability 




0.861 


0.048 


z 5* 3 


21(1) 


0.450 ±0.039 


0.167 ±0.034 


Z < 3 


274(13) 


0.244 ±0.008 


0.000 ±0.007 


Test Statistic 




5.905 


5.726 


Null-hypothesis Probability 




3.53 x 10-' 


1.03 x ur 8 



NOTE. — For the detailed definition of the Peto-Prentice test statistic see Feigelson & 
Nelson (1985). The null-hypothesis probability was calculated from the test statistic based 
on a Gaussian distribution, e.g., 1 — Pa = 7.27 x 10~ 7 , where Pa is the cumulative Gaussian 
probability at 4.954o\ 
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TABLE 7 

X-ray Properties of z > 4 HRLQs Derived from the Observed-Frame Ultrasoft Band (0.3-1.0 keV) 



Object Name Count Rate a fx,us b /2kcV,us c Oo x .us d Ag! ox ,rqq tjs c Aa OXjRL Q,us f a x.us - «ox g 

Chandra Cycle 12 Objects 



Archival X-ray 



SDSS J141209.96 + 062406.9 


1.32; 


-0.86 
-055 


0.81 


5.66 


-1.56 


0.13 


-0.11 


-0.05 


SDSS J142048.01 + 120545.9 


3 SQ+ 1 - 26 


2.48 


18.75 


-1.33 


0.34 


0.07 


0.01 


Data Objects 


















PSS 0121 +0347 


5 93+ 121 


2.02 


18.41 


-1.53 


0.22 


-0.02 


-0.06 


PMN J0324-2918 


6.62, „ 


4.37 


38.69 


-1.40 


0.35 


0.08 


0.00 


PMN J0525 - 3343 


34.33 


+0.87 
-0.85 


22.07 


178.82 


-1.15 


0.60 


0.33 


0.02 


Q 0906 + 6930 


5.32; 


-0.47 
0.43 


3.55 


31.09 


-1.26 


0.45 


0.18 


0.05 


RXJ1028.6-0844 


21.80 


+0.73 
-0.70 


16.25 


118.41 


-1.16 


0.56 


0.27 


-0.07 


CLASS J1325 + 1123 


2.01; 


-0.89 
-0.64 


0.75 


6.43 


-1.61 


0.11 


-0.13 


-0.08 


GB 1428 + 4217 


53.55 


+ 1.41 
-1.37 


32.45 


281.56 


-1.00 


0.73 


0.45 


-0.07 


GB 1508 + 5714 


22.33 


+0.51 
-0.50 


6.37 


48.61 


-1.03 


0.60 


0.27 


-0.07 



NOTE. — The objects included in this table are our z > 4 HRLQs with m.\ < 20, i.e., those z > 4 HRLQs used in the two-sample analyses in ^4. 1 1 except for the two 
Swift sources (PMN Jl 155-3107 and PMN J1951+0134). 

d The Chandra or XMM-Newton count rate in the observed-frame ultrasoft X-ray band (0.3- 1 .0 keV), in units of 10~ 3 s _1 . 

k The Galactic-absorption corrected flux in the observed-frame ultrasoft X-ray band, in units of 10~ 14 erg cm -2 s -1 . 

c The X-ray flux density at rest-frame 2 keV derived from the ultrasoft band count rate, in units of 10~ 32 erg cm -2 s -1 Hz -1 . 

^ The ctox values obtained from the ultrasoft X-ray band data. 

e The Aa ox rqq values obtained from the ultrasoft X-ray band data. 

The Adiffx blq values obtained from the ultrasoft X-ray band data. 
g The difference between the ce 0?t us values in this table and the a ox values in Column (15) of Table|4|which are obtained from the observed-frame soft band (0.5-2.0 keV) 
data. The X-ray spectral curvature discussed in j)4. 1.3 l is expected to result in negative values for this difference. 



TABLE 8 

Emission-Line REW Measurements 



Object Name MJD 11 REW(Lya + N V) REW(C IV) 

Chandra Cycle 12 Objects 



SDSS J102623.61 - 


-254259.5 b 


53734 


41.6 




SDSS J141209.96- 


- 062406.9 


53504 


2.4 


7.1 


SDSS J142048.01- 


- 120545.9 


53885 


55.0 


52.7 


SDSS J163950.52- 


-434003.6 


52051 


23.3 


7.8 


SDSS J165913.23- 


-210115.8 


52913 


17.9 


6.4 


[-ray Data Objects 










PSS 0121 + 0347 c 




53328 


81.1 


33.5 


SDSS J091316.55- 


- 591921. 6 b 


51907 


110.9 




SDSS J123503.03- 


-000331.7 


51941 


39.5 


54.1 


CLASS J1325 + 1123 


53148 


75.8 


63.4 


GB 1428 + 4217 d 




50283 


17.1 


12.9 


GB 1508 + 5714 




52079 


60.1 


47.0 


GB 1713 + 2148 d 




50280 


107.0 


127.7 



NOTE. — All REW values are in units of A. 
a The MJD (Modified Julian Date) listed here is the date of optical/UV spectroscopy. 
b Two objects (J0913+5919 and J1026+2542) do not have C IV coverage in their SDSS 
spectra because of their high redshifts. 

c Reference for emission-line REW measurements: Vignali et al. (2003a). 

d Reference for emission-line REW measurements: Hook & McMahon (1998). 
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TABLE 9 

Non-Thermal Dominance (NTD) Calculation 



Object Name 


logic iv 


logvL 13 50,obs 


log VLi350,pred 


NTD 




(erg s" 1 ) 


(ergs -1 ) 


(erg s _1 ) 




Chandra Cycle 12 Objects 










SDSS J141209.96 + 062406.9 


44.289 


46.670 


45.881 


6.150 


SDSS J142048.01 + 120545.9 


44.968 


46.465 


46.467 


0.995 


SDSS J163950.52 + 434003.6 


44.946 


47.255 


46.448 


6.402 


SDSS J165913. 23 + 210115. 8 


44.121 


46.516 


45.736 


6.021 


Archival X-ray Data Objects 










SDSS J123503.03 -000331.7 


44.597 


46.100 


46.148 


0.896 


CLASS J 1325 + 1123 


45.425 


46.860 


46.861 


0.997 


GB 1508 + 5714 


44.873 


46.503 


46.385 


1.312 



TABLE 10 
Correlation Analysis Results 





REW(Lya + N v) 




REW(C IV) 




N' 


Spearman 




Spearman 






rs l-fs 




rs 


AGIox.RQQ 
AOm.RLQ 
logi? 


32 
32 
298 


-0.19 71.0% 
-0.21 75.0%) 
0.14 98.0%) 


169 
169 
3536 


0.14 93.3% 
0.08 69.0% 
0.30 > 99.99% 



a N is the sample size. 



TABLE 11 

Joint X-ray Spectral Analysis Results 





Power Law 
with Galactic Absorption 


Power Law 
with Galactic and Intrinsic Absorption 


Object Name 


r x C/n a 


T x N H (10 22 cm- 2 ) C/n* 


J1026/J1412/J1420/J1659 b 
J0741/J1659 c 


1.59±S|5 122.47/148 
1.34+g;|| 46.87/66 


1.60±2-$ <9.50 122.47/148 
1.34+° ; *! < 12.41 46.87/66 



a C is the C-statistic, while n is the total number of spectral bins. 

b Our Chandra Cycle 12 targets with logR > 2.5. The names are in the format of 'ihhmm' for brevity. 

c Our Chandra Cycle 12 targets that are moderately radio loud. The names are in the format of 'Ihhmm' for 

brevity. 
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TABLE 12 

Upper Limits upon /-ray Flux and Luminosity from the Fermi LAT 



Object Name 



Photon Flux Limit 3 

1-9 



Energy Flux Limit 



(10 - * photons cm - s _1 ) (10 _1/ erg cm 



Flux Density Limit c 
(1(T 35 ergcm~ 2 s -1 Hz" 1 ) 



log(vti gcv) Limit d 
(erg s -1 Hz -1 ) 



Chandra Cycle 12 Objects 

SDSS J074154.71 + 252029.6 
SDSS J102623.61 + 254259.5 
SDSS J141209. 96 + 062406. 9 
SDSS J142048.01 + 120545.9 
SDSS J163950.52 + 434003.6 
SDSS J165913.23 + 2101 15.8 



4.5 
4.5 
5.0 
5.0 
4.5 
4.5 



2.4 
2.4 
2.6 
2.6 
2.4 
2.4 



2.1 
2.2 
2.0 
1.8 
1.6 
1.9 



47.4 
47.4 
47.3 
47.1 
47.1 
47.3 



Archival X-ray Data Objects 
PSS 0121 + 0347 
PMN J0324-2918 
PMN J0525 - 3343 
Q 0906 + 6930 

SDSS J091316. 55 + 591921. 6 

RXJ1028.6-0844 

PMN Jl 155 -3107 

SDSS J123503.03- 000331.7 

CLASS J1325 + 1123 

GB 1428 + 4217 

GB 1508 + 5714 

GB 1713 + 2148 

PMN J1951 +0134 



4.5 
4.5 
4.5 
3.5 
4.0 
5.0 
5.0 
5.0 
5.0 
4.0 
4.0 
5.0 
7.0 



2.4 
2.4 
2.4 
1.8 
2.1 
2.6 
2.6 
2.6 
2.6 
2.1 
2.1 
2.6 
3.7 



1.6 
1.9 
1.8 
1.8 
1.9 
1.9 
1.9 
2.1 
2.0 
1.7 
1.5 
1.8 
2.5 



47.1 

47.3 
47.2 
47.4 
47.3 
47.2 
47.2 
47.3 
47.3 
47.2 
47.1 
47.1 
47.3 



a The upper limit upon the Fermi LAT photon flux between observed-frame 100 MeV-100 GeV. 
k The upper limit upon the 7- ray energy flux between observed-frame 100 MeV-100 GeV. 
c The upper limit upon the y-ray flux density at rest-frame 1 GeV. 
The upper limit upon the monochromatic y-ray luminosity at rest-frame 1 GeV. 



TABLE 13 

Measurements of Black-Hole Mass and Eddington Ratio 



Object Name 



Chandra Cycle 12 Objects 

SDSS J141209. 96 + 062406. 9 
SDSS J142048.01 + 120545.9 
SDSS J163950.52 + 434003.6 

Archival X-ray Data Objects 

SDSS J123503.03 - 00033 1 .7 
CLASS J1325+ 1123 
GB 1508 + 5714 



log M BH (M e ) log (L/L Edd ) 



9.859+1.218 -0.708 
9.284 + 0.180 -0.338 
10.628 + 0.168 -0.893 



9.182 + 0.727 -0.602 
9.455 + 0.071 -0.115 
8.517 + 0.388 0.467 
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FIG. 1 . — The distribution of our sample of RLQs in the plane of M, (SDSS absolute i'-band magnitude) vs. logft (the logarithm of radio loudness), compared 
to high-redshift, moderately radio-loud quasars and RQQs. The filled stars show our Chandra Cycle 12 targets, while the filled squares are HRLQs at z > 4 with 
sensitive archival X-ray coverage. The asterisks show the high-redshift, moderately radio-loud quasars reported in Bassett et al. (2004), Lopez et al. (2006), 
and Mil. The small dots represent the high-redshift, radio-quiet SDSS quasars that have sensitive X-ray coverage. All the symbols are color-coded based on 
their redshifts using the color bar at the top right corner of the figure. The dotted and dashed lines show our criteria for RLQs (logR > 1) and highly radio-loud 
quasars (\ogR > 2.5). The dash-dotted curve shows the quasar radio-loudness distribution from Ivezic et al. (2004), which illustrates that our sample represents 
the quasars residing in the tail of high radio loudness. 
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FIG. 2. — The monochromatic luminosity at rest-frame 5 GHz (radio; upper panel) and at rest-frame 2500 A (UV; lower panel), plotted against redshift. The 
filled stars show our Chandra Cycle 12 targets, while the filled squares are our HRLQs with sensitive archival X-ray coverage. The plus signs represent the 
radio-loud and radio-intermediate objects in the full sample of Ml 1. The upper (lower) panel is color-coded based on UV (radio) luminosity using the color bars 
at the bottom right corner of each panel. Our sample of quasars are among the most luminous objects in both the radio and UV bands. 
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J1639 + 4340 

z = 3.980 

RQQ=-0.01 

log R= 1.621 




J0913+5919 
z = 5.122 
E-A« w RQQ =-0.09 
og R = 2.561 



J1235-0003 
:- z = 4.673 

4".,, rqq= 0-39 
og R = 2.928 



Lymari limit Lya 



C IV 





J1412 + 0624 
z = 4.467 

Aa o*, rqq= 0- 
og R = 2.698 





J1026+2542 
z = 5.304 

Aa ox, bqq= P- : 
log R = 3.542 




Lymari limit Lya C IV 

800 1000 1200 1400 1600 1800 
Rest-Frame X ( A ) 



600 800 1000 1200 1400 1600 1800 
Rest-Frame X ( A ) 

FIG. 3. — The SDSS spectra of five of our Chandra Cycle 12 targets and four objects with sensitive archival X-ray coverage, ordered by log/f which is shown 
at the top-left corner of each panel. Also labeled are the object names and redshifts. The y-coordinates are the flux density (Fx) in linear units. The spectra have 
been smoothed using a 20-pixel sliding-box filter. The wavelengths corresponding to major emission lines (Lya A 1216 and C IV A 1549) and the Lyman limit are 
labeled by the dotted vertical lines. The spectral resolution is R fa 2000. The apparent broad features of J1639+4340 and J1412+0624 at ~ 1000 A are caused 
by the stretched scaling due to their weak Lya emission lines. 
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FIG. 4. — The relation between radio loudness and a ox (top panel), Ao^ x R^rtmHdle panel), and AOo X rlq (bottom panel) for RLQs at all redshifts. The 
filled stars show our Chandra Cycle 12 targets, while the filled squares are our HRLQs with sensitive archival X-ray coverage. The open circles (downward 
arrows) represent the radio-loud and radio-intermediate objects in the full sample of Ml 1 that have X-ray detections (upper limits). The typical error bars for a ox , 
Agh ox ,rqq, and Aa„ x rlq for our objects are shown at the bottom-right comer of each panel. The dashed lines represent the positions of Aa ox = 0. The thick black 
lines show the mean Aa ox rqq values for the Mil objects binned in logi? (AlogS = 0.2 per bin). All symbols are color-coded based on their redshifts using the 
color bar at the top of the figure. 
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FIG. 5. — The histograms of Aa ox rqq (top panel) and Aa ox rlq (bottom panel) for the full-sample objects in Ml 1 with log/f > 2.5, z < 4, and m, < 20 (grey 
and open histograms for X-ray detected and undetected objects, respectively) and our HRLQs at z > 4 with m, < 20 (black histogram). 
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FIG. 6. — The REWs of the broad emission lines Lyce + N V (upper panel) and C IV (lower panel) as a function of radio loudness. In the upper panel, the open 
circles represent the radio-loud SDSS DR5 quasars in Table 1 of Diamond-Stanic et al. (2009) which gives the REW measurements of Lya + N V. In the lower 
panel, the small dots represent the radio-loud SDSS DR7 quasars with C IV REW measurements from Shen et al. (2011). Our objects with available line REW 
measurements are shown in both panels as filled stars (for Chandra Cycle 12 targets) or filled squares (archival objects). The weak-line quasar in our HRLQ 
sample (J1412+0624) is labeled in the upper panel. The typical error bars for emission-line REWs of our objects are shown at the bottom-right corner of each 
panel. All symbols are color-coded based on their redshifts using the color bar at the top of the figure. 
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FIG. 7. — The REWs of the broad emission lines hya + N V (top row) and C IV (bottom row) as a function of Aa ox rqq (left column) and Ag^rlq (right 
column), respectively. The open circles in the top-row panels represent the Mil objects with available REW(Lya + N V) measurements in Table 1 of Diamond- 
Stanic et al. (2009). In the bottom-row panels, the open circles are the Mil objects with REW(C IV) measurements from Shen et al. (2011). Our objects with 
available line REW measurements are shown in both panels as filled stars (for Chandra Cycle 12 targets) or filled squares (archival objects). The typical error 
bars for emission-line REWs, Aa ox .rqq, and Aa ox rlq of our objects are shown at the bottom-right corner of each panel. All symbols are color-coded based on 
their redshifts using the color bar at the top of the figure. 
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FIG. 8. — The relation between X-ray power-law photon index Tx and radio loudness log/?. The rilled stars show our Chandra Cycle 12 targets, while the 
rilled squares show our HRLQs with sensitive archival X-ray coverage. The open circles represent z > 4 moderately radio-loud quasars in Saez et al. (2011). 
The grey horizontal lines (rightward arrow) represent the 1 a error bars (lower limit) for the X-ray power-law photon indices. The dashed line shows our best-fit 
correlation between Fx and logi? for z > 4 RLQs, while the dotted line and dash-dotted line represent the Tx-logS correlations for z > 2 RLQs in Saez et al. 
(2011) and for z < 2 RLQs in Reeves & Turner (2000), respectively. The two moderately radio-loud objects among our Chandra Cycle 12 targets are apparent 
outliers (J0741+2520 and J1639+4340; see the two filled stars labeled in the bottom left part of the figure). 
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FIG. 9. — The rest-frame broad-band SEDs for our sample of objects ordered by logS (in descending order), including radio (red filled upside-down triangles), 
mid-infrared (WISE, orange open diamonds), near-infrared (2MASS and UKIRT, green open triangles), optical (blue filled circles), X-ray (purple asterisks and 
lines), and y-ray (filled square) data points. Downward arrows indicate upper limits. The purple solid (dotted) lines show the X-ray power-law spectra (and their 
uncertainty range) based on the photon-index values provided in Column 9 of Table [4] The purple asterisks represent observed-frame 2 keV. The black solid 
lines show the composite SEDs for the 10 RLQs in Shang et al. (2011; Sll) with comparable optical luminosity and radio loudness, scaled to the flux density 
at rest-frame 2500 A (corresponding to 10 15 1 Hz; see i]4.4t . The vertical dotted lines show the frequencies of rest-frame 2500 A and 2 keV. The name, redshift, 
and log R of each object are labeled at the bottom right of each panel. 
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FIG. [9] — Continued. 
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FIG. [9] — Continued. 



